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Abstract 

We consider measure-valued processes X t that solve the following martingale prob- 
lem: For a given initial measure X$, and for all smooth, compactly supported test 
functions tp, 

X t {<p) = X Q {<p) + \ [ X s (A^)ds + d f X s (ip)ds- f X s (L s <p)ds + M t (cp). 
^ Jo Jo Jo 

Here L 8 (x) is the local time density process associated with X t , and M t (tp) is a martin- 
gale with quadratic variation [M(<p)]t = J X s ((p 2 ) ds. Such processes arise as scaling 
limits of SIR epidemic models. We show that there exist critical values c (d) £ (0,oo) 
for dimensions d — 2,3 such that if 9 > 9 c (d), then the solution survives forever with 
positive probability, but if 9 < 9 c (d), then the solution dies out in finite time with 
probability 1 . For d — 1 we prove that the solution dies out almost surely for all values 
of 9. We also show that in dimensions d = 2, 3 the process dies out locally almost surely 
for any value of 9, that is, for any compact set K, the process X t (K) = eventually. 

Keywords: Spatial epidemic, Dawson- Watanabe process, Phase transition, Local ex- 
tinction. 
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1 Main results 

In a spatial 577? -epidemic, susceptible individuals in a population become infected at rates 
proportional to the current numbers of infected individuals in their vicinities. An infected 
individual remains infected for a certain period of time, and then recovers, after which he/she 
is immune to further infection. This paper concerns a stochastic model for such a spatial 
epidemic in which the infected set evolves as a finite measure-valued process Xt and the 
recovered set as the occupation measure Rt = J*q X s ds of the process Xt . If /i is a measure 
on M. d and ip is an integrable or non-negative measurable function on M. d , we will write n{ip) 
or (fj,, p) for the integral J pdfi.The dynamics of the model are specified by the following 
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martingale problem: for an initial value Xq = fi and any test function <p G C^(R ), 



M<p) 



^) + \ f x s (A<p)ds + e fx, 

z Jo Jo 



(ip)ds 



X s (L a <p) ds + M t {tp). (1.1) 



Here Mt(<p) is a continuous martingale with quadratic variation [jW(y)]t = J X s {ip 2 ) ds, and 
Lt(x) is the local time density process of Xt, that is, for each t > the function Lt(x) is the 
density of the occupation measure Rf. (Throughout this article, unless otherwise specified, 
the martingale Mt(<p) in a martingale problem such as (11.1)) will be a martingale relative to 
the minimal right continuous filtration of the process X, i.e., (J 7 ^) : = (r\ u> tcr(X s , s < u))). 
Dawson's Girsanov formula (Section 12. II below) implies that on a suitable probability space 
there exists a solution to (II. ID . that solutions are unique in law, and that the law is absolutely 
continuous with respect to the law of super-Brownian motion (on finite time intervals). 
Henceforth, we shall call a measure-valued processes Xt satisfying (II. ID a spatial epidemic 
process with transmission rate 6 and initial mass distribution ll. 

Measure- valued processes governed by the ma r tingal e problem (jl.ip arise as weak limits 
for discrete mode l s of S IR epidemics: see lLalley tod ) for the one dimensional case, and 



Lallev and Zhengl toi& i for the two and three dimensional cases. iMueller and Tribe! (|201ll 



suggest (see the remark at the end of Section 6 therein) that they sh ould also occur as 
scalin g limits of certain stochastic reaction diffusion systems. Earlier, Mueller and Tribe] 
(119941 ) studied processes governed by a related martingale problem in dimension d = 1, 
where the local time density in the third integral of (jl.ip is replaced by the density of Xt 
(this exists only in the one-di mensional case) . Su c h mo dels arose as the scaling limits of long- 
range contact processes (see IMueller and Tribe! (j 19951 )) which differ from the discrete SIR 
epidemics in that there is no recovered state where the individual is immune from infection. 
Unlike solutions to (jl.ip . the limiting models for the contact process were Markov processes. 

A measure- valued process Xt survives iff Xt(l) > for all t > 0. We say X dies 
out or becomes extinct iff Xt = for large enough t. For all of the above processes, the 
question of wheth e r or n ot there is survival or extinction is of funda mental importance- 
Mueller and Tribe! (120111 ) conjecture, by analogy with their results in IMueller and Tribe 
(1994J), that there is a critical value 6 C = 6 c {d) G (0, 00) for the transmission rate below 
which extinction is certain and above which survival has positive probability Our main 
result states that under a mild restriction on the initial measure \x this is true in dimensions 
d = 2 and d = 3, but that in d = 1 extinction is certain at all values of the parameter 9. 
The restriction on the initial measure is as follows: 

Assumption 1.1. The measure fj, has compact support and finite total mass, and when 
d = 2 or 3, its convolution fi* qt with the integrated Gauss kernel 



<lt(x) 



e -\x\ 2 /2t 

^p s (x)ds, where p t (x) = /2 , 



(1.2) 



is jointly continuous in t and x. 



Theorem 2 of ISudtanil (jl98fll ) asserts that in dimensions d = 2 and d = 3 a super- 
Brownian motion with finite initial mass distribution /x satisfying Assumption 11.11 has a 
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local time density process Lt(x) that is jointly continuous in t > and x 6 M . Since the 
law of a spatial epidemic is absolutely continuous relative to that of super-Brownian 
motion, spatial epidemics must also have jointly continuous local time processes in d = 2,3. 
In dimension d = 1, the existence and continuity of the local time process follows from the 
fact that the state of a super-Brownian motion at any time t is absolutely continuous with 
respect to Lebesgue measure, with a jointly continuous density. Thus, the equation (jl.ip 
makes sense in all dimensions d < 3, and so henceforth we assume that d < 3. 

Theorem 1.2. There exist critical values 9 C = 6 C (2) > and 9 C = 9 C (3) > such that 
the following is true: If d = 2 or d = 3 and Xt is a spatial epidemic process in M. d with 
transmission rate 9 and initial mass distribution /j, satisfying Assumption \l . 1[ then 

(a) if 9 < 8 C then Xt dies out almost surely, but 

(b) if 9 > 9 C then Xt survives with positive probability. 

If Xt is a spatial epidemic in M 1 with any transmission rate 9 and any finite initial mass 
distribution \i then X t dies out almost surely. 

Thus, in dimensions 2 and 3 a spatial epidemic can survive if the transmission rate is 
sufficiently high. The term — X s (L s ip) ds represents suppression of the disease as it tries 
to infect immune members of the population who have already been infected. This suggests 
that for the epidemic to survive it should move on to areas where the population has not 
yet been infected, that is, it should survive in transient waves which sweep through the 
population. In particular it should not survive locally as we now show is the case for any 
transmission rate 9. 

Theorem 1.3. Let Xt be a spatial epidemic with arbitrary transmission rate 9 6 K and 
initial mass distribution satisfying Assumption \1.1[ For any compact set K C K^, with 
probability one, 

X t {K) = eventually. (1.3) 

Consequently, with probability one the local time Lt{x) at any point x is eventually 
constant. Since the local time Lt{x) is jointly continuous in its arguments, it follows that 
-^oo^) := li m t-3.oo Lt(x) is finite and continuous in x almost surely. 

The proofs of Theorems 1 1 . 2Ifl~3l are somewhat intricate and technically involved, largely 
because of difficulties in carrying out comparison arguments for measure-valued processes 
defined by stochastic partial differential equations in which the entire histories of the so- 
lutions (e.g., local time density) influence the coefficients. However, the ideas behind the 
results can be explained, at least roughly, in simple terms. Consider first the assertion of 
global extinction in one dimension. If the epidemic process Xt were to survive with positive 
probability, then on this event its total mass Xt(l) would diverge to oo, since otherwise 
the process would be presented with infinitely many opportunities to become extinct (see 
Lemma r2.14p . In additio n, by a l arge d eviations calculation on a dominating super-Brownian 



motion with drift 9 (see iPinskyl (ll995l )). there exists c < oo such that Supp(Xt) C [— ct, ct] d 
for all large t. Therefore, for d = 1, on the event of survival and for large t, the average 
value of Lt(-) must satisfy 

/ct ft 
L t {x)dx = {2ct)- l \\L t \\ 1 = {2cty l / X s {l)ds^oo. 
-ct JO 
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If (jl.ip were valid for the function ip = 1 (it is only assumed for compactly supported 
functions), then it would follow that for large t the drift term in (jl.ip for the total mass 
X t (l) would eventually turn (very) negative, making it impossible for X t (l) to remain 
positive. The formal proof in Section 14,31 makes this heuristic argument precise. 

A local variation of this argument (which is harder to justify rigorously - see Section [7]) 
explains the strong local extinction asserted in Theorem 11.31 We will show that in order for 
Xt(K) > to occur at indefinitely large times, for some ball B centered at the origin, it 
must be the case that Xt(3B) integrates to oo. This, however, would imply that the local 
time in 2B \ B would grow indefinitely, eventually making the drift in the equation (jl.ip for 
X t (B) negative. 

A different line of argument makes it at least plausible that in dimensions d > 2 the 
epidemic Xf might survive with positive probability when the transmission rate 9 is suffi- 
ciently large. If 9 is large, then equation (jl.ip implies that when the infection first enters 
a region K of space it will, at least for a while, grow at least as fast as a super-Brownian 
motion with a large constant drift. Thus, with high probability, the total mass Xt(K) will 
become large long before the local time Lf becomes appreciable in K. In particular, for a 
cube K, if the size Xt(K) of the infected set reaches a certain threshold before the local 
time exceeds a fraction of this level, then the epidemic will have high probability of spread- 
ing to neighboring cubes quickly, and the infection in these cubes will have similarly high 
probability of spreading to neighboring cubes, and so on. Since high-density oriented site 
percolation in dimensions d > 2 has infinite clusters, with positive probability, it should 
then follow that the epidemic will reach infinitely many cubes with positive probability. It 
will take some work to implement this plan. This is done in Section [5] after some important 
groundwork is laid in Sections [2] and [3l 

For the extinction pa rt of Theorem 11.21 we will adapt the corresponding argument of 



Mueller and Tribe! (|1994l ) to our setting. For small positive 9 it is possible to rescale X so 



that the total mass process can be dominated by a subcritical branching process which dies 
out. The actual implementation of this idea is carried out in a slightly different manner in 
Section [6] (see Proposition 16. ID . A key observation, used here and elsewhere in this work, is 
that if the initial state is split up into pieces, then one can couple the epidemics so that the 
survival probability is dominated by the sum of the survival probabilities corresponding to 
the pieces (see Lemma |2. 18|) . 

The remainder of the paper is devoted to the proofs of Theorems ll.2til.31 The plan is as 
follows. In Section [2] we discuss existence and uniqueness of solutions to a class of martingale 
problems including (jl.ip . weak convergence of certain discrete processes to spatial epidemics, 
and basic comparison principles. In Section [3] we discuss some regularity properties of 
(supercritical) super-Brownian motions and their local time densities. In Sections 14.21 and 
4.31 we prove that the critical values 9 C in dimensions 2 and 3 do not depend on the initial 
mass distributions, and we prove that spatial epidemics in M 1 die out almost surely at all 
values of the transmission rate 9. In Section [5] we prove that spatial epidemics in dimensions 
2 and 3 can survive if the transmission rate 9 is sufficiently high; and in Section [6] we prove 
that at low values of 9 extinction is certain. We prove a weak form of local extinction in 
Section 14.11 and finally, in Section [7J we prove Theorem 11.31 

Our results extend to a somewhat larger class of measure- valued processes that incorpo- 
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rate location-dependent local suppression. Let K : M. d — > M + be a bounded, continuous (or, 
more generally, piecewise continuous), nonnegative function; call this the suppression rate. 
A spatial epidemic with local suppression rate K, transmission rate 8, branching rate 7 > 0, 
and inhibition parameter (3 > is a solution to the martingale problem (MP)^ 7 

X t (ip) = n{tp) + [ X s {A<p/2 + 0<p-K<p- PL s {X)<p) ds + yfyM t (<p) (1.4) 
J 

where Mt(<p) is a continuous martingale with quadratic variation [M(c^)]t = J Q * X s {tp 2 ) ds 
and Lt(x) is the local time density of Xp When K = and (3 = 8 = and 7 = 1, a process 
X satisfying (II. 4p is a super- Brownian motion: more generally, when K = and /3 = 0, 
it is a super- Brownian motion with drift 8 and branching rate 7; and when f3 = it is a 
super- Brownian motion with killing, with local killing rate K, drift and branching rate 7. 
Theorems 11.21 and 11.31 extend to all processes governed by (jl.4p with j3 > 0; the critical 
values # c will depend on the parameters, but not on the suppression rate function K if we 
restrict K to be compactly supported. In the interests of simplicity we shall prove our main 
results only in the case K = 0. However, solutions of the martingale problem fjl .4[) will be 
needed in the proofs of the main results even in the special case K = as they will arise 
naturally in the Markov property of solutions to (11.11) . 

Standing Notation. For any Borel subset D C W 1 , let M(D) be the space of finite 
Borel measures on D, equipped with the weak topology, and let Ai c (D) be the subset 
consisting of all measures with compact support in D. These spaces are partially ordered 
in a natural way: we write fx < v to mean that for all nonnegative, bounded functions cp, 

J Lpd[i < J ipdv. 

For a measure /x 6 Ai(D) and a nonnegative measurable function f : D —> M+, we shall 
continue to use the shorthand notations n{f) or (fj,, f) to denote the integral of / against \i 
and also write for /u(l), the total mass of fi. Finally, we let Cp^^R-f) be the space 
of non- negative, piecewise continuous, compactly supported and bounded functions on W d . 
(A function / will be called piecewise continuous if there exist finitely many bounded, 
simply connected domains £7j, each with piecewise smooth boundary, such that on each Ui 
the function / is bounded and continuous, and such that / = off Ucl(£/j), where cl(iTj) 
stands for the closure of Ui. A sequence /„ £ C p (M d ,M+) converges to / if and only 
if (a) the sequence f n converges to / uniformly on compact subsets of each Ui and (b) 
lim I ||in||oo — ll/llool = 0. The implied topology is metrizable, and the resulting metric space 
is separable and complete.) 

Conventions. Throughout the paper, C, C\, etc. denote generic constants whose values 
may change from line to line. Also, for any a, b £ R, a Ab := min(a, b) and a V6 := max(a, b). 
Finally, we use a "local scoping rule" for notation: Any notation introduced in a proof is 
local to the proof, unless otherwise indicated. 
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2 Preliminaries on the epidemic processes 



2.1 Dawson's Girsanov Theorem; existence and uniqueness 

Existence and uniqueness of soluti ons (in the weak se n se) to a class of martingale problems 



similar to fjl -4|) was established in iMueller and Tribd (|201ll ) using Daws on's G ir sanov the- 
orem. Existence i n the special case K = 0, 9 = was also proved in iLallevI (120091 ) and 



Lallev and Zhend (l2O10h . by weak convergence methods, which extend trivially to the gen- 



eral case. Nevertheless, since Dawson's Girsanov formula will be of crucial importance in 
many of the arguments to follow, we begin by reviewing the essentia l facts . We first state a 
variant of Dawson's Girsanov Theorem (Theorem IV. 1.6 in IPerkind (12002^ tailored to our 
needs. 

Let ft = D([0, oo); A4 c (R d )) be the canonical path space for compactly supported measure- 



valued processes, with coordinate maps Xt : Q — > M. C {R ) and associated filtration F = 
(J r ^)t>o- Fix a probability measure P on (f^J^), and suppose that there is a linear 
mapping tp i— > (Mt(ip))t>o from the space C 2 (R d ) to the space of F— adapted, continuous 
martingales such that Mo(ip) = and such that M(tp) has quadratic variation [M(ip)]t = 
f* (X„ f 2 )ds. This mapping extends to an orthogonal martingale measure dM(s, x) (see 
Walshl (fl98&k For any previsible x Borel process B : R + x $7 x R d — > R, we say that B is 



L 2 — admissible if 

[ (X s , B 2 S ) ds < oo for all * > P - almost surely. (2.1) 
Jo 

If B is L 2 — admissible, then the stochastic integrals 

r-t f 

B s (x)dM(s,x) (2.2) 



o 



exist and constitute a continuous, F— adapted local martingale with quadratic variation 
process Jq(X s , B 2 ) ds. Consequently, for each 7 > 0, the process 



Sf =exp (-L f I B M dM(s,x) - 1- f\x s ,B 2 ) ds) 
V\/7Jo JR d 2 7 Jo ) 



(2.3) 



is a continuous local martingale. 



Lemma 2.1. (Dawson's Girsanov Theorem) Let P be a probability measure on (£l,J-oo) such 
that under P the coordinate process (Xt)t>o satisfies the following martingale problem: For 
some uniformly bounded L 2 — admissible integrand A, 



ft ft 

Xtty) = Xoty) + £ / (X s , AV) ds + / (X s , A s ^) ds + ^M t p (ip), (2.4) 
z Jo Jo 



where Mf(ip) is a continuous -martingale with quadratic variation [M p (ip)]t = Jq X s (ip 2 ) ds. 

(a) Suppose that Q is another probability measure on (Q, J-oo) such that under Q the coor- 
dinate process Xt satisfies the martingale problem 

ft ft 

Xtty) = Xoty) + £ / (X s ,Aip) ds+ (X s , {A a + B s )^) ds + VlM?(iP), (2.5) 
z Jo Jo 
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for all if) 6 C 2 (M. d ), where B is a uniformly bounded L 2 — admissible integrand and M\ (ip) 
is a continuous martingale (under Q) with quadratic variation [M^(^)]t = X s (il) 2 ) ds. 
Suppose also that the restrictions of P and Q to the a— algebra J 7 ^ are equal. Then for each 
t < oo the measures P and Q on T% are mutually absolutely continuous, with likelihood ratio 



dQ 
dP 



=£t B - (2.6) 



In particular, Q is uniquely determined by the martingale problem (|2.5p . 

(b) Conversely, if Q is the probability measure determined by the likelihood ratios (|2.6p . then 
under Q the process Xt satisfies the martingale problem (|2,5p . 



We next apply the above to prove that the martingale problem (ll.4p is well-posed. We 
introduce the following notation: for each x € M. d , let Q r (x) be the closed cube of side 
length r centered at x, and, for notational ease, Q(x) := Q\{x). For any continuous path 
X t valued in M c (R d ), define 

ru v \ t (v \ tX( x r foXs(Q e (x))ds 
L(t,X,x) = L t (X,x) = Lf(x) = hmsup — ^ . 



When there is no confusion, we shall suppress the dependence on X and abbreviate L t (X, x) 
as Lt(x). If Xt is an adapted process on the filtered space (Q, F) then L(t,X,x) is nonnega- 
tive, nondecreasing in t, and V x B d — measurable, where B d is the Borel u-field on M. d and 
V is the previsible cr-field. 

Theorem 2.2. Let [i G M c (R d ) satisfy Assumption [7771 and let K £ C p (R d , R + ). For any 
fixed 9 £ I and 7 > 7 denote by P^ = -P^g' 7 the law of a super- Brownian motion with 
initial mass distribution [i, drift 6, and branching rate 7. 



(a) If X solves the martingale problem (II. 4ft with initial value Xq = \x then the law P^^k 

Hit 

Ef where B(s, u, x) = ~{K{x) + (3L{s, X, x)) (2.7) 



P^' 7 of X on the canonical path space is unique and given by 



dP, 



dP^ 



and dM(s,x) is the orthogonal martingale measure under P„. Conversely, if P^,k is the 
probability measure specified by (|2,7p then under P^^k the coordinate process Xt satisfies the 
martingale problem 



(b) The mapping (fJ,,K) 1— > P^,k is jointly measurable with respect to the appropriate Borel 
fields. 

(c) Under P^^k the local time process Lt(x) is jointly continuous in (t,x) and almost surely 
is the density of the occupation measure Rt = f^Xgds. 

(d) Under the measure P^k the process (X, L) is strong Markov, that is, for any J-^ -stopping 
time t, 

Pii,K(X T+ . £ C\F T ) = Px T ,K+f3L T {C) almost surely on {r < 00} for all C £ Too- 



8 



(e) For any pair K,K' of suppression rate functions, the probability measures Pa,K an d 
Pu K' are mutually absolutely continuous on J 7 ^ , with likelihood ratio 



dP, 



H,K' 



dP, 



expj-^^ J [K(x) - K'{x))dM(s,x) 

{K(xf - K' {xf)X s (dx)ds\ . (2.8) 



1 

"27 



Remark 2.3. Assertion (b) guarantees that if Xq and Kq are random and J-q— measurable 
then the random probability measure Px ,K i s J~o~ measurable. Similarly, if X T and K T 
are Jv— measurable then Px T ,K T is Jv— measurable. Moreover, since Px ,K * s a regular 
conditional distribution on the canonical path space given J-q, it follows from (d) that 
the strong Markov property holds when the initial condition Xq and the suppression rate 
function Kq are random. 

Remark 2.4. Since the local time density Lt is not uniformly bounded on finite time intervals, 
the exponential process Sf is not a priori a martingale. Part of the assertion of the theorem 
is that in fact £f is a martingale, and hence that (|2.7p defines a probability measure on T£ . 



Proof of Theorem \2. 21 (a) First we claim that any solution Xt to martingale problem (II. 4ft 
has the property that its local time density Lf(x,to) is bounded in (t, x) for t in finite 
intervals and for every t < 00 has compact support in x for almost every to. This follows 
because on some probability space a version of Xt can be coupled with a super-Brownian 
motion Xt with drift 9 and branching rate 7 such that Xt > Xt for all t > almost surely . 
See Proposition 12. 12l below. or alternatively the proof of Proposition IV. 1.4 in Perkins! ( 20021 ). 



Since super-Brownian motion has a continuous local time process with compact support in 
any finite time interval, by Sugitani's theorem, it follows that the process Xt also has a local 
time density L x (t,x) with the advertised properties. 

Unfortunately, we cannot directly apply the previous lemma to conclude (I2.7p . because 
hf is not uniformly bounded in to. To circumvent this problem we use a localization 
argument. Fix < b < 00, and consider the exponential process £j^ T ^, where 

r(6) = wf{t : nwx|J3 t (a;)| > b}. 

X 

By Lemma [2 .11 the process £^ T ^) 1S a mar tingale, and so under the probability measure Q b 
specified by equation (I2.7P (with the stopped exponential martingale as the likelihood ratio) 
the process Xt satisfies the martingale problem (jl.4p with K(x) + f3L(s,X,x) replaced by 
its stopped value. But the preceding paragraph implies that for each t, Q b {r(b) < t} — > 
as b —7- 00, and so for each t > 0, 

hm l{r(6) < t}) = => E^Sf) = 1. 

It follows that Ef is a martingale under P^ and that under the probability measure defined 
by (12 .7p the process Xt satisfies the martingale problem (II. 4p . 

(b, c) These are easy consequences of (a), the continuity of fx — > Pa, and Sugitani's theorem. 
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(d) It suffices to consider a finite- valued r. By (c), the local time Lj is the occupation 
density of Xt under P^k, s ° it follows that 

L T+ t(X, x) = L T (X,x) + Lt(X T+ ., x) for all (t,x) almost surely. 

If Q(u) is a regular conditional probability for X T+ . given T T , then it follows easily from this 
that almost surely under Q(u) the coordinate process satisfies the martingale problem (jl,4p 
with K replaced by K + L T . Therefore, by the uniqueness in law of solutions. Q(cj) = 



P 



) almost surely. The strong Markov property now follows. 



(e) This follows immediately from (a). 



□ 



Remark 2.5. Lemma [2.11 ho lds equally well on the larger space of continuous A4(M. d ) -valued 
paths (as in iPerkinsI (j2002l )) . The proof of Theorem 12.21 also holds on this larger space if 
you start with compactly supported initial conditions. That is, the solutions necessarily 
have compact supports for all t by the domination in (a). This slightly strengthens the 
uniqueness part and may be used implicitly below without further comment. The main 
reason for restricting to compactly supported measures is the use of Proposition 12.91 below 
in the proof of our main result Theorem 11.21 



2.2 Discrete epidemic models 

Measure-valued processes that satisfy the martingale problem (|1.1|) arise naturally as weak 
limits of discrete, finite-population stochas t ic mo dels o f spatial epidemics . Here we describe 
one such class of models, following Lalley ( 20091 ) and Lalley and Zheng ( 2010l ). Several of 
the couplings we shall develop later in the paper involving measure- valued epidemics will be 
constructed by first building corresponding couplings for discrete epidemics, then using the 
weak convergence of the discrete to the measure- valued processes to prove that they extend 
to the measure-valued setting. 

The discrete SIR-d epidemic models take place in populations of size N located at each of 
the sites of the integer lattice Z d . We shall call N the village size. Each of the N individuals 
(or particles) at a site x 6 Z rf may at any time be either susceptible, infected, recovered, 
or removed. Infected individuals remain infected for one unit of time, and then recover, 
after which they are immune to further infection. The rules governing the transmission 
of infection are as follows: at each time i = 1,2,... , for each pair (i x ,s y ) of an infected 
individual located at x and a susceptible individual at y, i x infects 
Pn{x, y), where 



s y with probability 



PN(x;y) =p 9 N {x;y) 



i + e/N a 

(2d+l)N 



if \y 



< 1 and = otherwise, 



(2.9) 



where \z\ is the Euclidean norm of z and 



a = a(d) = 2/(6 - d) 

is the critical exponent, see Theorem 1 in lLallev! (|2009h and Theorem 2 in lLallev and Zheng] 
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(j2010h . For the SIR-d model with village size N, define 



x ' - 1 • — set of infected particles at x at time i; X i (x) := (x)\; 

set of removed particles at x (at time 0); K N (x) := \K N (x)\; 

X^ (x) = # recovered particles at x at time n; (x) := (x)|; and 



Xf(x) 
K N {x) 



i<n 

Xf : = UajXf (x), I W := U x K N {x), and := UjR^(x). 



The Standard Construction. We now describe a way to construct this p r ocess using 
a percolation stru c ture. (This is a modification of the constructions in Lallev ( 20091 ) and 



Lallev and Zhengl (j2010l U We shall call this the Standard Construction. The percolation 



structure is a random graph with vertex set Z rf x {1, 2, . . . , N}; the vertex (x, i) represents 
the ith individual (or particle) in the "village" V x situated at location x G For each pair 
(x,i) and (y,j) of vertices whose spatial locations differ by at most 1 (that is, \x — y\ < 1) 
a pat— coin toss determines whether or not there is an edge between (x,i) and (y,j). (As is 
often the case in such constructions, it is useful, for comparison purposes, to assume that 
these coin tosses are realized using independent Uniform-[0, 1] random variables.) Thus, the 
resulting random graph Q = Q N has edges only between vertices in the same or neighboring 
villages. 

The spatial epidemic is defined by a deterministic algorithm on the random graph Q. 
Since the village size N is fixed in this algorithm, we shall omit all superscripts N in 
the specification of the algorithm. The colors green, blue, red, and black will be used to 
denote susceptible, infected, recovered, and removed vertices in each generation. For the 
0th generation, designate K(x) vertices at location x as black; the set of black vertices will 
not change during the course of the epidemic. Next, color Xq(x) vertices in V x blue, and all 
remaining vertices green. (Thus, in generation there are no red vertices.) Now define a 
time evolution as follows. In generation n + 1, the set X n+ i(x) of blue vertices will consist 
of all vertices that were green in generation n and were connected by edges of the random 
graph to blue vertices (that is, vertices in X n = UyX n (y)). Finally, all vertices that were blue 
in generation n become red in generation n + 1, and remain red in all subsequent generations 
(that is, R„+i(x) = E n (x) U X n (x)). 

The virtue of this construction is that all quantities of interest can easily be described 
in terms of the geometry of the random graph Q 1 = Q \ U y lK(y) obtained by deleting all 
black vertices from Q. The set X n (x) consists of all vertices at site x at distance n in the 
graph Q' from the set of vertices that were colored blue in generation 0. Similarly, the set 
M ra (x) consists of all vertices at site x at distance < n from the blue vertices in generation 0. 
The set U :r R 00 (x) of vertices that are ever infected during the course of the epidemic is the 
union of the connected clusters of the blue vertices of generation in Q' . It is immediately 
obvious from this that the recovered sets M n (x) are non-increasing in the initial condition 
K(-), and non-decreasing in Xo(-) and the transmission parameter 9. 

□ 

Denote by 

P" = (P n (x, y)) XtyeZ a = (P n (y - x)) w>y& * (2.10) 
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the transition probability kernel of the simple random walk on Z d , that is, P n = P*F n - 1 is 
the nth convolution power of the one-step transition probability kernel given by 

P 1 (sc,y) = l/(2d + l) for \y-x\<l, (2.11) 

and let a 2 = 2/ (2d + 1) be the variance of the distribution Pi(0, •). Let G n (x,y) be the 
associated Green's function 

G n (x) :=J2 P i( x )> 

and for any finite measure [i on % d denote by \iG n {x) = fi * G n (x) the convolution of \i 
with G n . Let (/j,Gt)(y) be the continuous extension to [0,oo) x M. d by linear interpolation. 

Next, we explain the re-scaling of the discrete epidemics that gives weak convergence to 
the measure- valued epidemics determined by the martingale problem (ll.4p . 



Definition 2.6. The Feller- Watanabe scaling operator Tn scales mass by 1/N a and space 
by l/VN a a 2 , that is, for any finite Borel measure fj,(dx) on W d and any test function (f(x), 



(<p, Fnim) = N~ a / ip(VN<*a 2 x) fx(dx). (2.12) 



Definition 2.7. The Sugitani scaling operator Sn scales mass by 1/N a( - 2 d l 2 ^> and space by 
1/V N a a 2 , that is, for any function K, 



K(\/N a a 2 x) 



The follo wing weak convergence theorem is a slight variant of the main results in lLallev 



(j2n09h and lLallev and ZheneJ (j2nid l. Let C b (R d ) (C c (R d ), respectively) denote the space 



of bounded continuous functions (compactly supported continuous functions, respectively) 
on M. d . When d = 1, we pose the following stronger condition on the initial measure [i: 

Assumption 2.8. The measure /i satisfies Assumption (jl.ip ; moreover, when d = 1, it has 
a density Xq(x) G C c (M). 

Proposition 2.9. Assume that d < 3, and suppose that the initial configurations are finitely 
supported measures fx N := Xq and K N such that for some measure fi satisfying Assump- 
tion \2.8\ either of the following conditions are satisfied: 

(a) if d = 1, then 



XtfiVN^ 2 x) 



'N c 



X {x) G C C (R); (2.14) 



(b) if d = 2 or 3, then 



F N ^ N => fi, (2.15) 
S N (M N *G [Naii )=>iJL*q t GC 6 ([0,oo) xR d ), (2.16) 

where the second convergence is in D([0, oo); Cb(M. )); and 
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(c) in all dimensions, 

S N K N => K G C p (R d ,R+). (2.17) 
Then we have the following weak convergence 

(TnX^c^SnR^^) ==>- (X t , L t (x)) (2.18) 

in D([0, oo); .M c (R rf )) x D([0, oo); C^R^)), where the limit process X has initial configuration 
Xq = [X, solves ([1.4p with 7 = 1, (3 = 1 and suppression rate K , and L t {x) is its local time 
density process. 

Remark 2.10. Proposition 12.91 implies, among other things, that for any initial mass dis- 
tribution fi G A4 c (M. d ) satisfying Assumption 12.81 and any suppression rate function K G 
C p (R d ,R+), the measure- valued epidemic process Xt satisfying (|1.4p is a weak limit of ap- 
propriately scaled discrete SIR epidemics. For given fj,, the initial conditions fj, N for the 
approximating processes can be constructed in the following manner. When d = 1, for each 
x G Z/VN a a 2 , let fi N = [y/N* ■ X (x)}. When d = 2 or 3, suppose that Supp(/x) C Q L (0) 
(the cube of side length L centered at 0), then for each N, divide Q ^ ^0,^2 ^(0) into non- 
overlapping cubes Qi of side length [N@], where a/2 > /3 > a(l/2 — 1/d), and within each 
cube Qi distribute := \p(Q%/V N a cr 2 ) x A^ a ] particles on the [N^] d integer sites in the 
cube in a uniform manner. One can then show that {fi N } satisfies (12. 14ft when d = 1 or 
(12.151) and (|2,16p when d = 2 or 3. The initial conditions K N can be constructed in a similar 
manner. 



Remark 2.11. The arguments of lLalley (j2009h and lLallev and Zhenel d20inh are based on 



the fact that each of the discrete SIR epidemics has law absolutely continuous with re- 
spect to the law of a critical branching random walk with the same initial condition. The 
Radon-Nikodym derivatives can be written explicitly as products, and these can be shown to 
converge to exponentials of the form Ef appearing in (|2.7p . Since branching random walks, 
after rescaling, converge to super-Brownian motions, it follows that the rescaled discrete 
SIR epidemics converge to processes related to super-Brownian motion by (|2.7p . that is, 
processes that solve the martingale problem (11. lj) . 

Routine modifications of these arguments can be used to establish weak convergence for 
a variety of discrete processes similar to or related to the discrete SIR epidemics constructed 
above. In particular, the convergence (I2.18P can be extended to joint weak convergence for 
coupled SIR epidemics with suitable initial conditions. For example, let [jl n,a , fi N,B be initial 
conditions satisfying the hypotheses fj2.15|) (12.16|) . and let ~& N,A , \ N ' B , be discrete SIR 
epidemics all constructed using the same percolation structure £7 , with the same initially 
removed sets ~K N (x), in such a way that the sets ~Kq' A (x) and ~Kq' B (x) are non-overlapping, 
with cardinalities fi N,A {x) and /j, N ' B (x), and such that 

X$(x) = X*> A (x) U X*' B (x). (2.19) 

Then after rescaling, the processes X N ' A , X N ' B , and X N converge jointly in law to (depen- 
dent) measure-valued epidemics X A ,X B , and Xt, with initial mass distributions fi A ,/j, B , 
and [i A + [i B , respectively, whose local time densities satisfy 

L A \I Lf < L t < L A + Lf. (2.20) 
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(The arguments that follow will not rely in an essential way on this joint convergence, 
however. All that is needed is that subsequences converge jointly, as this is enough to guar- 
antee the existence of coupled measure-valued processes satisfying the same monotonicity 
properties (such as (|2.20p ). Joint convergence along subsequences follows trivially from the 
weak convergence of marginals, since this implies tightness.) 

□ 



The discrete SIR-d process is naturally associated with a branching envelope. This 

— N — N AT 

is a nearest-neighbor branching random walk X n with initial condition X = X$ and 
offspring distribution Bin((2d + l)iV, p%(0, 0)) that dominates X^ , that is, such that for 
each n > and x 6 Z d , 

X*{x)<X N n {x). 



See Section 1.6 of iLallev! (l2009t ) for details concerning the construction. Since branching 



random walks, after Feller- Watanabe rescaling, converge weakly to super-Brownian motions, 
the vector-valued processes (X N ,X ), similarly rescaled, have marginals that converge 
weakly. It follows that after rescaling the laws of the vector- valued processes {X N ,X 1 are 
tight. Hence, any subsequence has a weakly convergent subsequence, and the limit process 
(X, X) must satisfy Xt < Xt- The component processes Xt and Xt of any such weak limit 
must be a measure- valued epidemic (that is, a solution of the martingale problem (|1.4|) with 
7 = 1) and a super-Brownian motion with drift 9, respectively. This implies the following. 



Proposition 2.12. Let Xt be a solution of the martingale problem (|l,4j) . Then on some 
probability space a version of Xt can be coupled with a dominating super-Brownian motion 
Xt > Xt, which we will call the super-Brownian motion envelope, with the same initial mass 
distribution fj,, and drift 9. 

Remark 2.13. Of course this has also been asserted in the proof of Theorem 12.21 However, 
the argument here provides an independent - and perhaps more natural - proof. 



2.3 Comparison lemmas 

The construction of the measure- valued spatial epidemic process as the weak limit of discrete 
epidemics and the Girsanov formulas (I2.7j) - (]2.8p lead to a number of basic comparison 
principles that will be used in the proof of Theorem 11.21 We formulate these as couplings, 
in which two epidemic processes (or super-Brownian motions) are constructed on a common 
probability space in such a way that various functionals of the processes (e.g., the limiting 
local time densities L OQ (x)) are ordered. 

Lemma 2.14. Suppose that Xt has law P^'jf 1 for some K <E Cp(M d ,M + ) and some initial 
condition fi that satisfies Assumption \2.8l Then 

P(X survives) = P I lim / Lt(x) dx = oo I = P( lim \XA = oo). 

Proof. This uses the existence of a coupling between the measure-valued epidemic Xt and 
its super-Brownian motion envelope. If Z s := \X S \ is the total mass at time s, then 
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linif^oo J* Kd Lt(x) dx = J °° Z s ds. Because Z s is continuous and is an absorbing state 
(e.g. by the strong Markov property in Theorem 12. 2p . if J Q Z s ds = oo then Xt must sur- 
vive. On the other hand, if limsup^^ Z t < oo, then there exists M £ N and an infinite 
sequence of stopping times r n — > oo such that 

T n +i > r n + 1, and Z Tn < M. (2.21) 

Consider the time period [r n ,r n + 1]. By the strong Markov property and the existence of 
a monotone coupling between a spatial epidemic and its super-Brownian motion envelope, 
the process Zt+ Tn is dominated by Feller diffusion with drift 9 and initial total mass less 
than M . This dies out in the next one unit of time with positive probability, independent 
of n, hence so does Xt+ Tn for t < 1. It follows that with probability 1, if X survives then 
lim^oo Zt = oo. As the latter trivially implies J °° Z s ds = oo, the proof is complete. □ 

Lemma 2.15. Fix 9 < 9* and 7 > 0. For any initial mass distribution fi £ _M c (R rf ) that 
satisfies Assumption \2.8\ and any K £ C p (M. d , there exist on some probability space 

epidemic processes (X,X*) £ D([0, 00); A4 c (M. d )) 2 with laws P^'jf 1 and P^^" 1 and local 
time densities Lt,L5, respectively, such that almost surely, for every t>0, 

L t < L* t . (2.22) 

Proof. This follows from the weak convergence result (|2. 18j) and the Standard Construction. 
Recall that in the standard construction of the discrete SIR epidemics, the evolution is 
determined by the random graph Q' in which edges are present with probabilities pn(x\ y) = 
p%{x;y) given by (|2.9p . These probabilities are increasing in 9. Consequently, it is possible 
(using auxiliary Uniform-[0,1] random variables) to simultaneously construct random graphs 
Q' and with percolation probabilities p e N (x;y) and p%(x;y), respectively, in such a way 
that the edge set of Q' is contained in that of Q\. This forces 

R n (x) C R* n (x) for all n > and x £ Z d , 

and hence also R n (x) < i?*(x). This inequality will be preserved upon taking weak limits, 
so we obtain (T2T22|> . □ 

Remark 2.16. It follows immediately from (|2.22p and Lemma [2.141 that 

P(X survives) < P(X* survives). (2.23) 

Lemma 2.17. Let K, K* £ C p (M d ;M + ) be suppression rate functions such that K < K* . 
Then for any fi £ A4 c (JH d ) satisfying Assumption ([2.8p . there exist spatial epidemics (X, X*) £ 
D([0, 00); A4 c (M. d )) 2 with marginal laws Pjjp and P^'][l and local time densities Lt,L\, re- 
spectively, such that almost surely, 

L t >L* for edit >0, (2.24) 
and P(X survives) > P(X* survives). 
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Proof. The existence of the coupling follows directly from the weak convergence (|2.18p and 
the Standard Construction, because in this construction, increasing the removed sets K(jc) 
decreases the sizes of the connected components. The assertion about survival probabilities 
follows from (|2.24p . by Lemma 12.141 □ 

Lemma 2.18. Let fj,,fj,Q,VQ E M c (M. d ) be initial mass distributions satisfying Assump- 
tion \2.8\ such that fj, = hq + vq. Then on some probability space there exist epidemic processes 
(Xt,Xl,Xf) E D([0, oo); Ai c (JH d )) 3 , with initial conditions fi,fiQ, and uq and marginal laws 
P q' 7 , Pu'oO ' > an d Pp'oQ respectively, such that 

uiax(L l t (x),L 2 t {x)) < L t (x) < L\(x) + L 2 t (x), for all t>0 and x G R d . (2.25) 
Consequently, 

P(X survives) < P(X 1 survives) + P(X 2 survives) and (2.26) 

P(X survives) > P(X l survives), for each % = 1,2. (2.27) 

Proof. This follows by the same argument as the preceding lemma (see Remark 12. lip . □ 

Lemma 12.181 describes the effect of adding infected mass at time 0. The next lemma 
concerns the effect of introducing additional infected mass at a time > after the epidemic 
has already begun. Let [i = fj,o + v$ be initial mass distributions satisfying the hypotheses 
of Lemma 12.181 Say that X% is a measure-valued epidemic with immigration at time i* if it 
satisfies the following martingale problem: For all ip G C 2 (R d ), 

Xt(<P) = M<f) + 1 [i,,oo)W^o(^) 

+ / X* (Ap/2 + 0<p - Kcp- PL*{X*)<p) ds + y/jM t *{cp), (2.28) 
J o 

where M 4 * is a continuous martingale with quadratic variation [M*((p)]t = Jq X*(ip 2 ) ds and 
L* = L x * is the local time density of X*. Existence and uniqueness of solutions to (|2.28p 
follows from Theorem 12.21 and the Markov property. 

Lemma 2.19. On some probability space there exist a solution Xt to (jl.4p with initial mass 
distribution = [1q + uq and a solution XI to (|2.28p such that 

Lf>Lf V t > and L* = . (2.29) 

Proof. This is by discrete approximation, using the Standard Construction of the discrete 
SIR epidemics. On each percolation structure Q = Q , we construct a pair of epidemics. 
The first, denoted by X = \ N , is constructed using initially infected sets Xq(x) = Hq (x) 
such that both (I2.15P and (|2.16p hold. The second, denoted by Y = Y , has initially infected 
sets Yq(x) C Xq(x) such that after Feller- Watanabe rescaling the initial mass distributions 
converge to ( see Remark 12.101 in Section 12. 2p . This second epidemic Y has spontaneous 
new infections at time [N a t*]: in particular, all individuals in the sets 

X \Y :=U X (X (x)\Y (x)) 
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who are not yet recovered become infected. Thus, the time evolution of the epidemic Y n 
is determined by the random graph Q' := Q\ UxK(x) as follows: (1) For n < [N a t*\, the 
recovered set consists of all vertices at graph distance < n from the initially infected 
set U2;Yo(x). (2) For n > [N a t*\, the set consists of all vertices v such that either the 
graph distance of v from L) x Yq(x) is < n, or the graph distance of v from Xo(x) \ Yo(x) is 
<n-[N a U}. 

By construction, for each N the set of individuals ever infected during the course of epi- 
demic X is identical with the set of individuals ever infected during the course of epidemic Y, 
because both are equal to the union of the connected components of all vertices in U y Xo(y) 
in the graph Q' . Moreover, in each generation n and at each site x the number of recovered 
individuals in the X— epidemic is no less than the number of recovered individuals in the 
Y— epidemic. Thus, the desired result will follow from the next assertion. Set 

Y t N =T N \Yf Nat] \ and 

where Tn and 5/v are the Feller- Watanabe and Sugitani rescaling operators. 

Claim 2.20. The vector-valued process (Y t N ,R^) converges weakly to a process (Xf,L^) 
such that X* solves the martingale problem (12. 28ft . and L* = L x is the local time density 
ofX*. 

Proof of the Claim (Sketch). This is done by fol lowing the likelihood ratio strategy described 
in Remark 12.111 As used in L alley ( 20091 ) and Lalley and Zheng ( 2010l ). this strategy was 



based on the fact that each discrete SIR epidemic considered had law absolutely continuous 
with respect to the law of a critical branching random walk with the same initial condition. 
The bulk of the proof consisted of showing that the likelihood ratios converged in law, under 
the branching random walk measure, to the Radon-Nikodym derivative of a measure-valued 
epidemic relative to the law of super-Brownian motion. For the processes considered in this 
claim, the appropriate comparison processes are not standard branching random walks, but 
rather branching random walks with immigration in which new particles are introduced at 
times [iV^i*] in such a way that after Feller- Watanabe rescaling the mass distributions of 
these new particles converge to uq. The laws of these processes converge, after rescaling, 
to the law of super-Brownian motion with immigration at time t*, that is, a process Y* 
satisfying the martingale problem (|2.28p with (3 = and K = 0. (This follows easily from 
the standard convergence theorem for critical branching random walks, because the effect 
of the immigration is simply to superimpose an independent branching random walk shifted 
in time by [N a U].) 

Consider the likelihood ratios for the law of the epidemic process Y t N relative to that of 
the corresponding branching random walk with immigratio n. Thes e are p roducts of factors 
indexed by (discrete) times t and lattice sites x £ Z d (see [Lalleyj ( 20091 ) . eqn. (53)). For 
t < [-/V a t*] the factors are exactly the same as in the case where there is no immigration. 
Beginning with time t = [N a t*\, new factors are introduced; these indicate the relative 
likelihood ratios for the newly introduced immigrants and their offspring. Under the law 
of the branching random walks with immigration the immigrants and their offspring evolve 
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independently of the progeny of the original (time 0) particles. Using th i s fact , one can 
show, in much the same manner as in lLallevI (l2Q09h and lLallev and Zhend (|2010h . that the 
likelihood ratios converge weakly (under the branching random walk with immigration laws) 
to the Radon-Nikodym derivative of the process X* relative to super-Brownian motion with 
immigration. In carrying out this final step, the main hurdle is showing that in the epidemics 
with immigration, the numbers of individuals in the sets Xo \ Yo who are infected prior to 
time [iV ^] is of order Op(l). Here is a brief synopsis of the argument: Since the local time 
densities, after rescaling, converge, the maximum number of recovered individuals at time 
[iV ^*] at any site is of order Op(N a ( 2 ~ d / 2 ^). Consequently, because Hq\Yq has cardinality 
on the order N a , the number of individuals in Xo \ Yo infected prior to time [iV ^*] is of 
order 

Op(iV«(2-rf/2)) 



0{N a ) 



N 



Op(1), 



since a = 2/(6 — d) . 



(2.30) 



□ 



This completes the proof of Lemma 12.191 □ 

In the proof of Theorem 11.21 it will be necessary to compare the evolution of a measure- 
valued epidemic X with a coupled process in which additional infected mass is introduced 
at a random time. Say that X^ is a measure- valued epidemic with immigration at time r if 
it satisfies the following martingale problem : For all ip 6 C 2 (M. d ), 

Xt(<P) = M<P) + 1 [r,oo)(i)^0(^) 

+ j X* s (A<p/2 + 9^ -Kip- pLf(X*)<p) ds + vWM, (2.31) 

where r is a stopping time relative to the filtration T x * , M t * is a continuous martingale 
with quadratic variation [M*(ip)]t = J * X*(ip 2 ) ds, and L is the local time density. 

Lemma 2.21. Let jj, = /io + vq be initial mass distributions satisfying Assumption \2.8[ 
Then on some probability space there exist epidemics (Xt,X^) G D([0, oo); Ai c (M. d )) 2 such 
that (1) Xt solves the martingale problem (II. 4ft with initial value Xq = n; (2) XI solves the 
martingale problem ([2.311) : and (3) 

L?>Lf V t > and L* = L*\ (2.32) 



Proof. By the usual continuity (weak convergence) arguments, it suffices to prove this for 
stopping times r that take values in a finite set. By a routine induction on the cardinality 
of this finite set, it suffices to consider stopping times that take values in a two-element 
set {0,t*}. For such stopping times, the result follows from Lemma 12.191 since this can be 
applied conditionally on Jo. □ 

Now the results of Lemmas 12.171 [2~T8l and 12.211 can be combined, allowing us to couple 
the measure-valued epidemic Xt with a measure-valued process X£ in which the infected 
mass is decreased and the suppression rate increased at a random time r. Here we will use 
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the strong Markov property (Theorem 12.2( c)) and Remark 12.31 The process X.\ will satisfy 
the following martingale problem: for every tp £ C 2 (M. d ), 

Ho(tp) + J X* (Atp/2 + 0<p-Kip- fiLf V) ds + ^M*{^), for all t < r, 

r T ( V ) + y" X s * (A^/2 + 9<p - K*<p - (3L X V) + VWfa) - M;^)), for all t > r, 

(2.33) 

where is the local time density and 

(i) r is a finite T x * -stopping time; 

(ii) Y T is an T x _ —measurable random measure satisfying Y T < X*_ + vq; 

(iii) K* is an J- x * —measurable random element of C p (K d ,K+) satisfying K* > K; and 

(iv) Mf(ip) is an T x * -continuous martingale with quadratic variation [M*(ip)]t = Jq X*(ip 2 ) ds. 

Proposition 2.22. Let fx = [xq + vq be initial mass distributions satisfying Assumption \2.8[ 
Then there exist epidemics (X t ,X^) G D([0, oo); M c (R d )) 2 such that 

(i) Xt solves the martingale problem (II. 4ft with initial value Xq = fx; 

(ii) XI satisfies the martingale problem (|2,33p : and 

(iii) the local time densities of X and X* satisfy 

L? > Lf* for all t > 0. 

This proposition can be used iteratively, using the strong Markov property and standard 
arguments, so as to allow immigration and increases in the suppression rate K* at increasing 
stopping times = To < T\ < t<i < ■ ■ ■ < oo. The associated martingale problem is as follows: 
for every ip 6 C 2 (W d ), and for all r, < t < 73+1, 

xiiv) = (h(<p) 

+ I X* s (A<p)/2 + 6<p-K*<p-p(L s (X*)-L Ti (X*))<p) ds + yfy(M?(<p) - M*^)), 
J n 

(2.34) 

where 

(i) for each i = 1, 2, . . ., fa and i>i are T x _ —measurable random measures such that 

fii + Vi< X*._ + 

(ii) Kq = 0, and K* £ C p (IR c( ,lR+) is, for each i, an T x —measurable random function 
such that 

Kt>Kt_ 1 + P{L n {X*)-L 7i _ l (X*)); and 

(iii) Ml is a continuous T x —martingale with quadratic variation [M* (tp)]t = Jq X* (<p 2 ) ds. 

The existence of a solution to this martingale problem follows from the strong Markov 
property (Theorem 12.21 (d)). Proposition 12.221 and a standard induction argument now yield 
the following comparison result. 
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Proposition 2.23. Let fj, = /io be initial mass distributions satisfying Assumption \2.8l 
Then on some probability space there exist measure-valued processes X and X* such that (i) 
X solves the martingale problem (MP) 'q specified in (|1.4j) ; (ii) X* solves the martingale 
problem (|2.34p .- and (Hi) the corresponding local time processes satisfy 

L? > L?* for all t > 0. 



2.4 The Sandwich Lemma 

Lemma 2.24. For any measure \x G Ai c (M. d ) satisfying Assumption (|2.8p . any k > 0, 

9 G M., and any function K G C p (M d ,M4-) there exist, on some probability space, measure- 
valued processes X,X,X_ with common initial state Xq = X = X_ = fi such that 

X t <X t < X t for all t < r, 

where 

t = inf |t > : max (K(x) + hf (x)) > k| , 

with the following laws: 

(i) Xt is a spatial epidemics with local suppression rate K, transmission rate 9, branching 

rate 7 = 1 and inhibition parameter (3 = 1; 
(ii) Xt is a super- Brownian motion with drift 9; and 
(Hi) X_ t is a super- Brownian motion with drift 9 — k. 

The proof will once again be based on discrete approximations. We shall build approx- 
imating discrete epidemic processes that satisfy the analogous sandwich relationship. The 
construction makes use of the following lemma. First, observe that in a discrete SIR epi- 
demic, when two infected individuals simultaneously attempt to infect the same susceptible 
individual, all but one of the attempts fail; call such an occurrence a collision. (Hence, 
for example, when three infected individuals simultaneously attempt to infect the same 
susceptible individual, then the number of collisions would be = 3.) 

By slightly modifying the proof of Lemma 9 in lLallev and Zhend (|2O10h . in particular, 



by noticing that the statement right above equation (62) therein also holds for the way that 
we count the number of collisions here, we get the following 



Lemm a 2.25 (A slight variant of Lemma 9 and equations (62) - (64) in lLalley and Zheng 



( 2010h ). For each pair (n, x) G N x let (x) be the number of collisions at site x and 



time n in the SIR-d epidemic with village size N. Assume that the hypotheses (I2,15!) -( |2.16p 
of Proposition \2.9\ are satisfied, for some a < 2/(6 — d). Then for any fixed T > 0, 

E E Y. T nW = °( Na )- (2-35) 

n<N a T x 

A direct consequence of the previous lemma is that if we define a Modified SIR process 
in the following way: 

Modified SIR process: At any site/time (x,t), each particle produces Bin((A r — K N (y) — 
(y)),p 9 N (x,y)) number of offspring at neighboring sites y, where (y) = ^2 i<k X^(y), 
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then the Modified SIR process can be constructed together with the original SIR process, 
in much the same way as for the branching envelope with the original SIR process, such 
that: (1) the Modified SIR process always dominates the original SIR process; and (2) the 
discrepancy Dt(x) > between them at site x and time t satisfies 

t ^ T ^ A(x)= ° p(iVQ) - (2 ' 36) 

X 

Therefore after the Feller- Watanabe scaling as in Proposition 12.91 the modified SIR process 
will converge to the same limit as in Proposition 12.91 

(To show flUSD , observe that if we let D n := Y, x D n (x), and D n = D n /(1 + 9/N a ) n , 
then D n is a sub-martingale: E(D n +i\J- n ) > D n , with the inequality due to collisions at 
generation n+1. Further note that for any T > 0, 

E ( D [N<*T]) < (l + 0/jV")[^ 

= °\ E \ E E^l } =0{N°). 



^n<N a T x 

(|2.36p then follows from the Doob's martingale inequality.) 
We now prove Lemma [2.241 

Proof of Lemma \2.24\ We shall build approximating particle systems that satisfy the analo- 
gous sandwich relationship. Choose Xq and K such that (j2.15j) ~ (j2.17p are satisfied. The 
super-solution X t (x) is a nearest-neighbor branching random walk with initial configura- 
tion Xq and such that at any site/time (x, t), each particle at site x produces Bin(N,p s N (x, y)) 

number of offspring at neighboring sites y. By Watanabe's Theorem, X^ (x) converges to 
the desired Xf{x). As noted above, the Modified SIR Xj^ will approximate Xp Define the 
stopping time 



t = min It > : max 



lax (K N (x) + R?(x)) > kN^-W} . 



The sub-solution Xf(x) is a nearest-neighbor branching random walk with initial config- 
uration Xq and such that at any site/time (x,t), each particle at x produces Bin([iV — 
KN a ( 2 ~ d / 2 ^],p e N (x,y)) number of offspring at neighboring sites y. It is clear that before 

time t n , the three processes xf (x), X^(x) and Xf (x) can be built on a common proba- 
bility space such that 

X? < X t N < Xt for all t < r N . 

By Watanabe's Theorem X N converges weakly to the super-Brownian motion X_ and taking 
weak limits in the above we complete the proof. □ 



21 



2.5 Scaling 



It will be necessary, in some of the arguments to follow, to rescale time and/or space. When 
a super-Brownian motion, or more generally a solution to the martingale problem (jl.4p is 
rescaled, its diffusion rate may change, that is, the Laplacian in (|1.4p may be multiplied by 
a constant a. The resulting martingale problem is as follows: 

X t (<p) = X (<p)+~ f X s {&ip)ds+0 f X s {ip)ds- [ (X s (K<p)+pX s (L s v))ds+y/jM t (<p), 
1 Jo Jo Jo 

(2.37) 

where a, (3, 7 > and 9 £ R are constants, K £ C p (M d ,R + ), and Mt(ip) is a continuous 
martingale with quadratic variation [M(c^)]i = X s (tp 2 ) ds. As usual, Lt is the local time 
density of the process Xf. We shall refer to this martingale problem as (MP)^ 7,a and 
continue to write (MP)^ 7 if a = 1. 

Lemma 2.26. Let X t solve the martingale problem (MP) 9 ^ l3 I f' a . For any constants a,b,c> 0, 
define a new measure-valued process Ut by 

Ut(tp) = c / ip(bx)X at (dx), for all bounded measurable ip on M d . (2.38) 

J X 



Ql Rl 



Then Ut solves the martingale problem (MP)" , K '? ,a with parameters 



r = ad, p' = ^£ : j' = acj, a' = ab 2 a, K'(x) = aK(x/b), 



and initial measure defined by f ip{x)n'{dx) = cf ip{bx) fi(dx). The local time densities 
L = L x and L u are related by 



L u s {x) = abd L as[-r) , for all x,t. (2.39) 
Proof. This is by routine calculations. □ 



3 Preliminaries on (supercritical) super-Brownian motions 

In this section we present some regularity results for super-Brownian motions and their local 
times. The results are only of interest, and only will be used, for d > 1 and so we assume 
d = 2 or 3 throughout this Section. 



3.1 Uniform regularity of super-Brownian motions 

In this and the following subsection, let Yj = be a (driftless) super-Brownian motion 
with initial state Yq = ^i, and let P^ = P^'q' 1 be its law. Denote by B r {y) the open Euclidean 
ball in M. d centered at y of radius r, and for any measure \x £ A4(K. d ), define 

D{», r) = sup{M(S r (y)) : y £ R d }. (3.1) 
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For any function cp and any measure [i G A4(R d ), set [lip = n*tp, where * denotes convolution. 
In particular, for any t > 0, 



(m)( x ) = Pt(x - y)fJ-(dy), and 0<7t)(x) 



p s (x - y)n{dy) ds, 

where pt and qt are the Gauss kernel and the integrated Gauss Kernel in (jl,2p respectively. 
For r G (0, 1] let 

/i(r) = ^ r log — , and (/?(r) = r 2 ^1 + log — ^ . 
Finally, for T,ro,Ci > introduce the event 

G T (r ; Ci, C 2 ) = {L>(y t , r) < d (D(fj,p t , C 2 r) + p(r)) for all r < r and t < T} . 



The following lemma is an easy consequence of the proof of Theorem 4.7 in lBarlow et al. 



(|199lh . 

Lemma 3.1. If K > 1 i/iere are constants C\,C 2 > (depending on K), and for any T > 
i/iere is an vq{K,T) 6 (0, 1] suc/i i/iai /or a// A > 1 and any ji with \fi\ = A, 



pjG T (roe~ x ; C^d)) >1 



-if A 



(3.2) 



Proof. This is a quantitative version of Theorem 4.7 of iBarlow et al.l (|l99ll ). The proof of 
that result shows for K > 1 there are constants Ci,C2,C3 > 1 such that for all A > 1, 
T > 0, n £ N and ^ with |//| = A, 



P^G T (h(2- n ); d,C 2 ) c ) < C 3 (T + 1)(A + 1)2- Kn + C 3 A2- An < C 3 (2T + 3)A2" An . (3.3) 

Here one has to c hase constant s a bit to check that the constant C2.2 in the proof of the 
above Theorem in I Barlow et al.l (|l99ll ) may be taken to be as large as you like at the cost 
of our C 2 and their C4.2 b eing large. The latte r can then be handled in the key bound in the 
proof of Theorem 4.7 in IBarlow et al.l (|199ll ) by taking our C\ large enough. Now choose 
no > 2 in N so that 



C 3 (2T + 3)2- Kno < e~ KX \- 1 < C 3 (2T + s)2~ Kno+K . 
The above definition implies 



(3.4) 



h(2 



2- n °/ 2 (n log2) 1 / 2 > 



2^C 3 (2T + 3) 



[21og2] 1/2 >r (K,T)e-\ 



where in the last inequality we used the fact that for all X,K > 1, A -1 /* 2 *) > A" 1 > 
exp(-A/2)/2. Therefore and imply that 



Pu{G T {r Q e-\C u C 2 y) < e 



-K\ 



□ 



23 



To formulate the next result we introduce the following 

Definition 3.2. For any positive constants A, A, and ro, with r$ < 1, and any measure 
/x G M(M. d ), we say that /i is (A, X,ro) — admissible if 

D(At,r) < Ar 2 (Ar d ~ 2 + (l + log = ip(r) for all r < r e~ A . 

Corollary 3.3. For any fixed K > 1 anc? T > ; f/iere exist positive constants A = A(K, T) 
and ro = ro(K,T) < 1, such that for all A > 1 and /i ureid |/i| = A, 

P M (Y T is (A, A, r )- admissible) > 1 - e" KA . (3.5) 

Proof. This follows from Lemma [3.11 by noticing that 

D(fi PT ,C 2 r) < C 3 (T)Xr d . 

□ 

3.2 Local time densities of super-Brownian motions 

Recall (equation (jl,2p ) that pt{x) and qt(x) are the Gauss kernel and the integrated Gauss 
kernel, respectively 

Lemma 3.4. Suppose that /i G A / l(M d ) satisfies = A and is (A, A,ro) — admissible for 
some constants A and ro. For any < /3 < 2 — d/2 and any fixed T > 0, define 

Qi(T) := max (/igr) (x), and Q 2 {T) := max / / s~^^ 2 p s (x — y) fJ.(dy) ds. (3.6) 
x x Jo J 

Then there exists a constant A' = A'(A,ro,T, f3) > such that for both i = 1,2 and for all 
A > 1, 

{ A'X 2 when d=2 
Qi(T) < K d (X) := < (3.7) 
I A'AV w/ien d = 3. 

Proof. We shall only prove the result for Qi(T); the proof for Q2(T) is similar. Let ro = 
roe _A . We first deal with the integral for t G [rj^ 3 ,T]: to do so, for any fixed x G M. d we 
cover M. d with balls B{ of radius ro and with center of distance di = kro to x for some 
k G Z>o- Then for any t > 0, 

2 ' 



(2^/ 2 ( m )(*) = |t-d/2 eX pf-^n ^(dy) 



< 



The balls can be chosen in such a way that for any k > 3 there are at most C(k — 2) d 1 
balls with center of distance kro to x. It is then easy to see that there exist constants Cj 
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such that 



^ /2 exp(-^^) n(dy) 

OO f 

<dt- d ^(7 ) + t~ d/2 exp ( 

1 — q \ 



(k - l)'rg 



2t 

k=3 



• (k - 2) d ~V(ro) 



/oo / 2' — 2 \ 



< 



CiV'^o)*" 1 + C 3 (A + (1 + log(l/f )) 2 ) when d = 2 

CiV(r )r 3/2 + C 3 (A + (1 + log(l/f )) 2 /^o) when d = 3. 



Therefore 



< 



C 4 ^(ro)(logT + log(l/f )) + C 3 T (A + (1 + log(l/f )) 2 ) when d = 2 

C 4 ^(f )f ( - 4/3) + C 3 T (A + (1 + log(l/f )) 2 /f ) when d = 3. 

which can be bounded by Kd(\) for all A > 1 for an appropriate choice of A' . Now we deal 
with the integral for t G [0, r Q ]: 



JO \J|a;-S/|<i 3 / 8 J\x-y\>t 3 / s J \ Zt J 

-8/3 ~S/3 



<C 5 + C7 6 A. 

□ 

The following lemma is implicit in ISugitanJ (Il989h . 

Lemma 3.5. Suppose that /j G .M(M d ) and /or a// i > 0, Qi(t) < 00, and /or some 
< /3 < 2 — d/2 and a// 1 > 0, Q2{t) < 00, where Q\ and Q2 are defined in (|3.6p . Define 

Z t (x) = L t (x) - (fj,q t )(x). 

Then for any T > 0, i/iere ezzsi constants 770 = Vo(fi,T) > 0,d = Ci(f3,T) > suc/i t/iai 
/or all < r] < t]q and t <T, 

£ M exp ( ^f^f W) ) < exp (dQ 2 (2i) • 77) , for all \a - b\ < 2; (3.8) 

and 

Ea exp (77 Z t (a)) < exp (C 2 Qi(2t) • n) , for all a G R d . (3.9) 
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Proof. The second claim (|3.9p follows from Lemma 3.4 in ISug itanil (jl989h . 

To prove (|3.8p . following (3.34) in lSug itanil dl989h . for a random variable X we say that 



Eex.p(rjX) = exp y^c w ??" 

\n=l 

holds formally if for all /c > 1, < oo and 

d k (exp ( En=l C" 7 ? 



r;=0 



By (3.38), (3.45) and (3.48) in ISugitanil (|l989h . we have formally 

^exp(7 ? (Z t (a) - Z t (b))) = exp (f)" (M, v n (t, ■))) , 

where for n > 2, and x S 

>>2t 

kn(t,x)| <6 n .| a -6|^t 2 -^)/ 2 / s-^ 2 (p s (a-x)+p s (b-x)) ds, 

Jo 

and {6 n } are defined inductively as follows: 

n-l 

bi = C A > 0; 6 n = C 5 ^ b k b n - k . 



(3.10) 



(3.11) 



fc=i 



Using the proof of Lemma 3.4 in ISugitanil (119891 ). if we let f(r]) = Yl^Li ^nV" then for some 
<5>0, 



f(r,)-C 4 r ] = C 5 f(rj) 2 ; f(r,) 



1-y/l- 4C 4 C 5 r? 
2C 5 



<Crj for < r] < 5. (3.12) 



This shows that Y^ n b n r\ n has a positive radius of convergence, and the formal equation (|3.10p 
is indeed an equation when r\ is sufficiently close to because the Taylor series for the analytic 
function on the right-hand side is given by the left-hand side. Relation (|3,8p then follows 
easily from the upper bounds (|3.1ip and (|3,12p . 

□ 

Corollary 3.6. Under the assumptions of the previous lemma, for any fixed T > 0, there 
exist constants r/o = t]q((3, T) > 0,Ci = Cj(/3, T) > such that for all < r] < t]q, 



E^exp 



( ^ |Lr |a-fe|fl T(6)l ) < 2 ^v{CiQ2(2T) n), forall\a-b\<2; (3.13) 



and 



exp (r/Ly(a)) < exp (C2<5i(2T) rj) , for all a G 



(3.14) 
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Proof. Relation (|3. 14[) follows easily from (|3.9p . As for (|3.13p . by fj3.8|) and the elementary 
inequality e\ x \ < e x + e~ x , 



( V \L T (a) -L r (6)| ^ < 2exp {CiQ2i2T)rj) . exp 



r/|(/i«r)(a) - (mtr){b)\ 



By (3.44) in Sugitani ( 19891 ) we have for all x, y and i > 0, 

\Pt(x) -Pt(v)\ < c(P)r^ 2 \x-yf{p2t{x)+ P2t {y)), 

and so 

livutrXa) - (jjutr)(b)\ <C\a-bf f [ t^l 2 {p 2t {a - x) + p 2t {b - x)) »(dx) ds 

JO Jx 

<C\a-bfQ 2 (2T). 

□ 

Lemma 3.7. Suppose that T(x) is an almost surely continuous random field on W 1 such 
that for some r] > and f3 > 0, 



Pexp ( ^ T f°Lfr|/?^ ) - Cl ' a// ' a ~ 6 ' - ^ and 

Pexp (r/T(a)) < C 2 , for all a £ R d . 

Let Ml = max ag Q i ( ) T(a). T/ien /or all L £~N and m > 0, 

,2d//3 i / - , ; / r ? 7Ti 



(3.15) 



P(M L > m) < (C ie za/P + C 2 )L a exp 



1+7 



(3.16) 



where 7 = 8d^^ 2 . In particular, for any K > 0, i/iere exists C > 0, depending only on K, 
C\, C 2 and j3, such that for all L > 1 anc? A > 1, 



P ( M £ > C ±±^L) < e~ K " 



(3.17) 



Proof. It suffices to prove f|3. 16|) for L = 1 as th e result s then follow trivially by divid- 
ing Ql(0) into unit cubes. We apply Lemma 1 of iGarsial ( 19721 ) with p(u) = \l/(it) = 

and 



exr, 1 JlM. 



dxdy. It is easy to check that this B satisfies the hypothesis 



of Lemma 1 in the above reference. That result, or more precisely (10) in the proof, implies 



< T(0) + - flog(P) + 2d / log(l/n) d(u p ) 
V 1 Jo 

= T(0) + ^[log(5) + (2d//3)]. 
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Therefore for x > 0, 

P(Mi > (1 + i)x) < C 2 exp(-r]x) + P(log(P) > r/x - (2d//3)) 
< C 2 exp(-r ? x) + £(J3) exp((2d//3) - t/s) 
= (C 2 + Cie M//3 )e- ??x , 

which is (|3.16|) for L = 1, and where (|3.15[) is used to see that E(B) <C\. □ 

Combining Lemma 13,41 Corollary 13.61 and Lemma 13.71 we obtain the following for the 
local time L of the super-Brownian motion Yj. 

Proposition 3.8. For any T > 7 M > 0, K > 0, A > anc? ro > there exists a constant 
A", depending only on (T, M,K,A,ro), so that for all A > 1 and all (A, A, ro) — admissible 
\i G satisfying = X, the local time, Lt{x), of the super-Brownian motion Y^ 1 

satisfies 

P„ ( max L T (x) > A"\ K d (X)) < e" x \ 

\\x\<Me x J 

where Kd(X) is defined in (|3.7p . 

Proof. By Lemma 13.41 and Corollary 13.61 for any fixed 77 < 770, if we let 77(A) = r]/kd{X), 
then for some fixed C\ and C2, for all A > 1, the assumptions (13. 15|) hold for the random 
field Lt(x) and < /3 < 2 — d/2 by replacing rj with 77(A). The conclusion then follows 
from (1317)1 . □ 

3.3 Local time densities of supercritical super-Brownian motions 

In this and the following subsection Yf = Yf" is a super-Brownian motion with drift one 
starting at an initial state /i, let P^ = P^'g' 1 be its law. Further denote by P^ = P^'q' 1 the 
law of a (driftless) super-Brownian motion starting at [i. By Lemma [2.11 we have 



:=$ t =exp (V°(l)-^Vsl^) 



where M % denote the martingale measure under P*, i = 0, 1, and therefore 

Af?(l)=M t 1 (l)+ / |Y,| da. (3.18) 
j 



Lemma 3.9. For any if > 1, T > 0, A > 1 and /j G A4(lR d ) wt/i = X, 

P° ($r • l ( $ T > e A-A) ) < — . (3.19) 
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Proof. Using (|3.18p we see the above expectation equals 



E° a $ T • 1 / \ I =Pl MUl) + - / \Y s \ds > KX 

i fr^i 



1 



^((Mf(l)) 2 ) | ^(Jo 



if 2 A 2 if A 

4A(e r - 1) A(e T - 1) 5e^ 
if 2 A 2 + KX ~ K 



□ 



Proposition 3.10. For any fixed T > and any e > 0, i/iere exist constants A = 
A(T,e) > and ro = ro(T,e) G (0, 1] suc/i that for all X > 1 and a// /i m'i/i |//| = A, 

P* (Y T is (A, X, r )- admissible) > 1 - e. (3.20) 

Proof. Let Gt,x denote the event in (|3.20p . For T,e as above choose K = K(e,T) > 2 so 
that 

^ + e -Kl2 < e 

and then choose A and ro as in Corollary 13, 3l for this choice of K and T, so that they depend 
ultimately on T and e. Then the previous Lemma and Corollary 13.31 imply that 

P 1 ,(G^ x ) = E^ t -1 G c tx ) 

< E°^ T ■ l ( * T > e A*/ 2) ) + e XK / 2 P%(G c T ,x) 
<We T /K + e- XK / 2 <e, 

where the choice of K is used in the last inequality. □ 

The same reasoning, but now using Proposition 13.81 in place of Corollary I3.3[ gives the 
following: 

Proposition 3.11. For any positive constants T, M and e, A > and ro > there 
exists a constant A", depending only on (T, M, K, A, ro), so that for all X > 1 and all 
(A, A,ro) — admissible jjL E Ai(M, d ) satisfying \fi\ = X, the local time, Lt(x), ofYj 1 satisfies 



Pi ( max L T (x) < A" X kJX)) > 1 - e. 

M \\x\<Me>> J 



3.4 Propagation of supercritical super-Brownian motions 

We continue to let Yt = Y^ be a super-Brownian motion with drift one starting at and let 
P^ denote its law. Recall that for x G Ij d , Q r (x) denotes the cube of side length r centered 
at x, and Q(x) := Qi(x). 
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Lemma 3.12. For any T > 1 and e > 0, there exists a constant M = M(T,e) > such 
that for any A > e and any /i satisfying = A and Supp(//) C Q(0), we have 



PMSupp(Y[0,T])CQ MVI ^(0) >l-s 



Proof. This is a direct consequence of Theorem A in iPinskvl (|1995l ) . □ 



Write W(0) = {x G Z, d : = 1} for the nearest neighbors of the origin in Z d . Fix a T 

sufficiently large such that 

min min e T (lofri * vt) (v) > 2, (3.21) 
where pt(x) is the Gauss kernel in ([1.2p . and * denotes convolution. 

Lemma 3.13. For any e > and T as above, there exists Ao = Xq(T,e) > such that for 
any \i satisfying Supp(/z) C Q(0) and = X > Ao, 

p l {Y T {Q{x)) > A for all x G Af(0)) > 1 - e. 



Proof. By a well-known moment formula (see, e.g., Exercise II. 5. 2 in IPerkind tooi )). to 



gether with the assumption that Supp(/x) C Q(0) and (|3,2ip . for any x £ A/"(0), 



EY T (Q(x)) = e T (/j, l Q(x) * PT ) > 2| M | = 2A, 



and 

rT 



Var(y r (Q(x))) < e 2T < M , / (l 
Consequently, by the Chebyshev inequality, 



o 



P(Y T (Q(x)) < A) <P (\Y T (Q(x)) - EY T {Q{x))\ > ^EY T (Q(x)) 

■ ^(Vjo i 1 Q(x)*P(T-s)) 2 *Ps ds) 



< 



(2A) 2 



<c T \-\ 

The conclusion follows. □ 

4 A weak form of local extinction and its consequences 
4.1 A weak form of local extinction 

Let Vd = 7r c( / 2 /r(l + d/2) be the volume of a unit d-dimensional ball. 
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Proposition 4.1. There exists k < oo such that for any 6 G R, 7 > 0, and i^T G (^(R^, R + ) ; 
i/X solves (MP)^ 1 ^ and satisfies Assumption ] then for any N > 1, 

E{Loo, 1 Bn (o)) < + V d (n + 29 + )(N + l) d , (4.1) 

and 

E(LI,1 Bn{0) ) < 4|//| + V d (K + 26+)\N + l) rf . (4.2) 

Proof. First, observe that there exists k > such that for any JV > 1 there exists a function 
if = (fx G C 2 such that 

\/\ip\ < and l Bjv (o) < V < iBjv+iW- ( 4 -3) 

For example, set <p{x) = ip(\x\) where if; = ipw '■ R —> [0,1] is a smooth, even function 
bounded above and below by the indicators of [— JV — 1,JV + 1] and [— N,N], monotone 
on [— JV — 1,—N] (and therefore also on [JV, JV + 1]), and such that (for instance) tp(x) = 
(x - N - l) 4 for x G [-JV - 1, -JV -1/2]. Because ip G C 2 (R d ), the martingale identity 
(jl.4p applies (with /3 = 1), so after taking expectations we obtain 

1 /** 

£(X t , 99) = (n, ip) + -E(L t ,A<p) + 9E(L t , <p) - E(L t , Kip) — E (X s , L s ip) ds. 

1 Jo 
A routine integration by parts shows that 

/"* 1 

J (X s ,L s ip) ds = -(L 2 ,ip). 

Since \Aip\ < and (Xt,ip) > 0, it follows that 

-2( / u, 99) < K,E(Lt, v^) + 20+£(L 4 , p> - £(L?, y>) (4.4) 

< (K + 29 + )E{L t ,^)-E{Llip) 

< (« + 2#+)£<L t , ^) - (KT^JV + l)- d ) (£<X t , ^)) 2 , 

the last by Cauchy- Schwartz and the fact that tp has support contained in Bn + i(0). This 
clearly gives an upper bound on E{Lt, y/ip) that is independent of t. In fact, completing the 
square on the right side and letting t — > 00 yields 

£<Loo, y^> < + V d (n + 2#+)(JV + l) d . (4.5) 

Relation (|4.ip follows, since yJTp bounds the indicator function of Bn(0). Finally, by the 
second inequality in (I4.4p . 

EiL^ip) <2\ l i\ + (K + 29+)E{L 00 ,y/Ip). 

Relation (0~2]) follows from (Oj) . □ 

Remark 4.2. The above Proposition easily shows that each of the terms on the left-hand 
side of (jl.ip converges a.s. as t ^ 00. Therefore Xf(ip) converges a.s. as t — > 00 and clearly 
the limit must be by the above. This shows that Xf(K) approaches as t — > 00 for all 
compact sets K a.s. Our Theorem 1.3 asserts a much stronger result, namely that Xt(K) = 
for large enough t a.s. 
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4.2 Universality of the critical values 9 C 



For any \x satisfying Assumption 12. 8| if X solves (jl.ip with 9 < 0, then P(X survives) = 
because X is dominated by a critical super-Brownian motion (by Proposition 12.12] ). which 
goes extinct almost surely. Lemma [2 .151 and Remark 12 . 1 61 therefore imply that for any such [i 
and any function K £ C p (M. d ; M_|_) there is a critical value 9 c (/j,, K) 6 [0, oo] so that a spatial 
epidemic X with suppression rate K and transmission parameter 9 (see (|l-4j) ) survives with 
positive probability if 9 > 9 c (fj,, K) and with zero probability if < 9 c (/j,, K). 



Proposition 4.3. The critical value 9 c (n,K) depends only on the dimension d and not on 



the choice of ^ ji satisfying Assumption ] 1.1\ or K € C p llc ''■ 



Proof. In this argument 9 will be fixed and 7 = 1, so we suppress the dependence of 
the laws P^]f on 9 and 7. By Theorem 12.21 for any measure \i 6 Ai c (R d ) satisfying 
Assumption 11.11 and any two suppression rate functions K,K' £ C p (R rf ; M_|_), the laws P^k 
and Pu,K' are mutually absolutely continuous on , with Radon-Nikodym derivative (|2.8|) . 
Since X and AT' both have compact support, inequality (|4.2p of Proposition 14.11 implies that 
the integrals in the likelihood ratio converge, so that 

K m ( d lt£L) :=Y 

t~>co \ dP^K ) Tt 

exists and is positive P^^k— almost surely. Hence, by Fatou's lemma, 

P^K'(X survives) = Hm P^ K ,{\X t \ > 0) > E^ K (Yl {x survives} 

It follows that if X survives with positive P^ t K probability then it also survives with P^^k' 
probability. Reversing the roles of K and K' shows that the reverse is also true. Therefore, 
the critical value 9(fj,,K) does not depend on K. 

To complete the proof, it suffices, in view of the preceding paragraph, to prove that if 
X survives with positive probability under P^q then it survives with positive probability 
under P V) q for v 7^ 0. By the Markov property (Theorem 12.21 (d)), 

P^oiX survives) = E^ q (P Xi ,l 1 (X survives)), 

and similarly for P v ,o- By the argument of the preceding paragraph, 

Pxi,L!(X survives) > <J=^ Px x fi{X survives) > 0, 

so for both ui = \i and u = u, 

Pui,o(X survives) > <J=^ E U! fi(Px 1 fi(X survives)) > 0. 

But the laws of X\ under P^o and P V) $ are mutually absol utely continuous. (This c an be 
seen as follows. First, by the absolute contin uity results in Evans and Perkins ( 199 ll ) (see, 



for example Theorem III.2.2 in&2 »), if P„ and P v are the laws of super^ownian 



A 4 

motions with initial conditions \i and u, then the distributions of X\ under P^ and P v are 
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mutually absolutely continuous. Second, by Theorem 12.21 fa), for any initial measure oj the 
measures P^ and P w ^ are mutually absolutely continuous.) Therefore, 

P fJi fl(X survives) > <J=^ P V $(X survives) > 0. 

□ 



Note that the above arguments do not require [i to satisfy the stronger Assumption 12.8 
instead just the original Assumption 11.11 



4.3 Extinction in dimension one 

Proposition 4.4. If d = 1 then for every 9 £ M and every initial measure fi that satisfies 
Assumption li.il the solution Xt of the martingale problem (jl.ip dies out almost surely. 

Proof. By the previous proposition, we only need to prove the conclusion for the /i's which 
satisfy the stronger Assumption 12.81 namely, those with compactly supported continuous 
densities. Fix any such an initial measure /i. Below we will also fix the mass creation 
rate 9, so we will use the abbreviated notation X = X^ ,d for the solution to the martingale 
problem (jl.ip . 

First, by Proposition 12.12) on some probability space there is a version of the process X 
and a super-Brownian motion Xt with drift 9 such that Xq = Xq = fi and Xt < Xt for all 
t > 0. 

By a result of Pinsky ( 19951 ). there is a positive constant C = Cq < oo such that almost 



surely the support of the random measure X is eventually contained in the interval [—Ct, Ct]. 
Since X dominates X, the same is true for X. Now by Lemma I2.14) on the event that X 
survives, the total mass of the measure Xt must diverge. Because this mass is (eventually) 
contained in [—Ct,Ct], it follows that the occupation density process Lt(x) must satisfy 



Ct 

L t (x) dx — > oo 

-ct 



on the event of survival. Hence, if there is positive probability of survival then 



i r l 

-E / L t (x) dx — > oo. 
t J~ct 



But this contradicts fj4. 1 j) in Proposition 14.11 □ 



5 Proof of survival when d = 2 or 3 

In this section we prove that in dimensions 2 and 3, for all sufficiently large values of the 
transmission rate 9, spatial epidemics - that is, solutions of the martingale problem (jl.ip 
- survive with positive probability. By Proposition 14. 3( the critical value 9 C for survival in 
dimensions d = 2, 3 does not not depend on the initial mass distribution ji, hence it suffices 
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to prove that for some finite measure [i there is positive probability of survival. The proof 
will make use of an auxiliary 3— dependent site ■percolation process: this will be constructed 
in such a way that if percolation occurs with positive probability then the epidemic must 
survive with positive probability. We will show that by taking 9 sufficiently large, we can 
make the density of the site percolation arbitrarily close to 1. Since percolation occurs with 
positive proba bility in a site percolation process when the density is near 1 (Chapter 4 of 
Purred (jl995h ). it will follow that for large values of 9 t he epidemic pro cess will survive with 



positive probability. We refer the reader to Chapter 4 of Durrett ( 19951 ) for terminology and 
a general framework for such comparison arguments. 



5.1 Scaled process 

We assume d = 2 or 3 throughout this section. Let Xt be a spatial epidemic process with 
transmission rate 9 and initial mass distribution fj,, that is, a solution to the martingale 
problem (jl.ip . It will be convenient to work with a rescaled version of the spatial epidemic 
defined as follows: for any 9 > 0, 

Utty) = 9X t/e (i>(V9 •)), for all e C*(R d ). 

The effect of this rescaling is described by Lemma [2.26l in particular, Ut satisfies the martin- 
gale problem (MP)l'J 1 with (3 = 6»( d " 6 )/ 2 and fl defined by J ijj(x) dfl(x) = 9 J ip(V9x) dfj,(x). 
For notational ease, we will use the notation 

P = p(0) = 9 {d ^ )/2 (5.1) 

in this section, and we will drop the tilde on the initial measure ji. We will show that when 
9 is sufficiently large, for a suitable initial condition /i, the process Ut survives with positive 
probability. 



5.2 Sandwich Lemma 

By Lemma 12.241 a spatial epidemic process can be bounded below and above by super- 
Brownian motions with different drift terms up to the time that its local time density 
exceeds some threshold. We now explain how the result of Lemma 12.241 translates to the 
rescaled processes. 

For any fj, G M. c (M. d ) satisfying Assumption 12.81 an d any function K G C p (M d ,IR + ), 
let U be a solution of the martingale problem (MP)^ 1 , i.e., the spatial epidemic with 
transmission rate 1, branching rate 1, inhibition parameter f3, local suppression rate K, and 
initial mass distribution [i. In addition, for any fixed constant k > 0, let Ut and U_ t be 
super-Brownian motions with drift 1 and drift 1 — /3 ■ k respectively. Denote by Mt and M_ t 
the orthogonal martingale measures associated with Ut and U_ t , respectively. 

Lemma 5.1. Versions of the processes Ut, Ut andU_ t , all with the same initial condition \i, 
can be built on a common probability space in such a way that 

U t <U t < U t for all t < r, 
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where 

t = inf |t > : max (K(x) + /3L t (U, x)) > /3k} . 

Proof. This follows by first rescaling U, U and U_ as in Lemma [2. 261 so that the (3 parameter 
becomes 1 and the drift parameters (or transmission rates) of U and U become 0(1 — (3k) 
and 9, respectively. Then one may apply Lemma 12.241 to the rescaled process. Finally 
undoing the scaling leads to the required conclusion. □ 

By Lemma [2.1 f a), the law of U_ is absolutely continuous with respect to that of U, and 
the likelihood ratio on Tt is 

LR t K = exp (-/3rcM t (l) - [ \U S \ ds\ . (5.2) 



5.3 Percolation probability estimates 

Recall from Section [2] that Q r (x) denotes the closed cube of side length r centered at x, 
and, as before, we abbreviate Q{x) := Q\(x). The auxiliary site percolation processes will 
be constructed by partitioning the space M. d into cubes Q(x) of side length 1 centered at 
lattice points x G Z d , and then using the behavior of the superprocesses in the cube Q(x) 
to determine whether the site x will be occupied or not in the auxiliary percolation process. 
Roughly, a site x will be occupied if, within a certain fixed amount of time T < oo, the 
measure- valued process U started from a certain initial mass distribution supported by Q(x) 
manages to generate a sufficiently large total mass in each of the adjacent cubes Q{y) while 
simultaneously not accumulating too much local time. The objective of this section is to 
develop estimates that will allow us to conclude that if 6 is large (hence f3 is small), then x 
is occupied with high probability. 

Define the grid V to be Z+ when d = 2 and Z+ x {0} when d = 3, where Z+ := {x = 
(x\,X2) G Z 2 : Xi > 0, i = 1,2}. For x,y G F, we say that 

| if ||d|i < ||y||i or 

[ IMIi = II2/II1 and x i < yi- 

Here ||x||i := \xi\ is the L 1 -norm. This defines a total order on T, and so the points of 
the lattice can be enumerated as = x(l) -< x(2) -<.... The notation x ^ y is understood 
as x -< y or x = y. Define A(x) to be the set of y E T such that x -< y and \\x — y\\i = 1. 
For example, when d = 2, A(0) = {(0, 1), (1,0)}. We shall call any y G A(x) an "immediate 
offspring" of x, and x an "immediate predecessor" of y. 

Fix T so large that (j3~2T|) holds. For any e > 0, let A(e) = A(T,e) and r (e) = r (T,e) 
be the constants specified in Proposition 13 . 1 Ol For any measure- valued process Xt with local 
time density Lt(X), and for any xGT, M > 0, x>0, A>0 and e > 0, define the following 



35 



events: 

f F\M;X,x) = {Supp(L r pO) C Q M {x)}; 

F 2 (x;X) = {m a xL T (X,y)< X }; 

v (5.4) 
F 3 (X,x) = {X T (Q{y)) > \X \, for all y G A(x)}; 

F 4 (e;X) = {X T is {A(e/4), \X \, r (e/4))-admissible}. 
Observe that these events depend on the choice of T. For brevity we will write 

F l {M) = F^M; 17,0), F 2 ( X ) = F 2 ( X ;U), F 3 = F 3 (77,0), and f\e) = F\e;U). 

Define functions 

/ (0) - { 6,1/2 When ^ = 2 (5 5) 

I log when cf = 3, 

and 

M = M(M,9) = My / logf d (9), and X = x(A",0) = A" f d (0)K d (f d (9)), (5.6) 
where K d (-) is the function defined in (|3.7p . 

Lemma 5.2. For any eo > , there exist positive constants 9q,M and A", depending only 
on T and eo, such that if > 6q then for any initial measure \i supported by Q(0), of total 
mass \fj,\ = fd(6) an d (A(£q / 4), f d (9),ro(eo/ 4:)) -admissible, the super- Brownian motion Ut 
with drift 1 and initial mass distribution fx satisfies 

p(F 1 (M)nF 2 ( X )nF 3 nF i (e )) >1-b . (5.7) 



Proof. This is a direct consequence of Lemma I3.12[ Proposition 13.111 Lemma 13.131 and 
Proposition 13.101 More specifically, by Lemma |3.12[ there exists constant M = M(T,£q) 
such that P[F (M)) > 1 — eo/A. Moreover, by Proposition 13.11] there exists A" such that 
P(F 2 (x)) > 1 - £o/4. (Note that in Proposition EOH A" only depends on (T,M,£,A,r ), 
and in our case the M, A and ro all only depend on (T, eo), so ultimately A" only depends 
on (T, £q).) Next, by Lemma 13.131 there exists 9q > such that for any 9 > 9q, 

P(F 3 ) > 1 - e /4. (5.8) 

Finally, by Proposition 13.101 as long as 9 is such that f d (9) > 1, P(F 4 (eo)) > 1 — £q/4. 

□ 

The next result explains the choice of f d {9). 

Corollary 5.3. For any positive constants 9,M and A", let U_ t = ]J^ K be the super- 
Brownian motion with drift 1 — (3k and initial mass distribution jj,, where 

k = [M + lfx = (Ms/\ogf d {9) + l) 2 • A'f d (9) Kd (f d (9)). (5.9) 

Then for any £q, there exists 9q > such that if 9 > 9$ and if the initial condition fi is 
supported by Q(0) and of total mass = f d {9), then 

P (F 3 (U, 0)) > 1 - 3e /2. (5.10) 
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Proof. Let F 5 (e ) = {LR^ < 1 + e } for the likelihood ratio LR£ defined in (|5.2j) . Using 
the fact that on the canonical path space (with coordinate maps Ut) Af*(l) = M t (1) — 
(3k fo \LLs\ds, we have 

E U (LR% • l( F 5( eo) )c) 

= pH (-/3kM t (1) - (f3 2 K 2 /2) £ \U s \ds > log(l + e )j 



J |E^|ds>log(l + e )V 



Here and below we use E u and P u (E— and P-, resp.) to indicate that the expectation 
and probability are taken with respect to the law of U (U, resp.). We have that 

t-T 



I \Hs\ ds < f d (0) [ e s ds< C T f d {e). 
Jo Jo 



This and the definitions of (3 (in (|5.1|) ). K and imply that E—((3 2 k 2 \U S \ ds) = o(l) as 
6 goes to infinity. It now follows easily from Markov's inequality and the above that there 
exists 6>o > such that for any 9 > 6q, E u {LR 1 ^ ■ l(F 5 (a )) c ) < £ o- Therefore by (|5.8p . the 
complement of the event F 3 (U_, 0) has probability bounded above by 

□ 



e + E u \ LRPf • <i+ £() } " 1 (F 3 y) - e ° + ( X + e °) ' T - •> 



Combining the sandwich lemma (Lemma 15. ip and the previous two results we obtain 

Proposition 5.4. For any eq > f/iene exisi positive constants 0q,M,A" such that, for 
any 9 > 0q, any initial condition /i satisfying the hypotheses of Lemma \5.2[ and any K £ 
C p (lR d ,R + ) such that 

K(x) + p X ■ Wo) (*) < /3k, /or a// x G M d , 



the process U solving i^MP) 1 '^ satisfies 



p(F 1 (M;U,0)nF 2 (x;U)nF 3 (U 1 0)nF\e ;U)) > 1 - 3e , (5.11) 



where the events F % (i = 1,2,3,4) are defined as in ([5.4p by replacing X with U and x 
with 0. 

Proof. On the event ^(MjnF^x), 

L T (U,x) < X -1q~{o)(x)- 
Therefore by the assumption on K and Lemma [5.11 



U t <U t < U t , for alii < T on F^M) n F 2 (x). 

The required bound now follows from Lemma [5. 2 [ Corollary 15 . 31 and an elementary argument. 

□ 
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5.4 Proof of survival 



Proposition 5.5. For some finite measure \i and some 6 < oo, if Ut solves (MP) q with 
j3 = #( 6 ~ d )/ 2 , then 

P(U survives) > 0. 

Proof. Fix a T so that (|3,2ip holds. Fix Eq > small enough such that any 3-dependent 
oriented site percolation process on Z+ with density at least (1 — 6eo) has positive probability 
of percolation. For this Sq, let 6 > 8q, where 8q is as in Proposition 15.41 Then choose a 
measure \x so that it satisfies the hypotheses of Lemma 15.21 with £q specified as above. Let 
Lt(x) denote the local time density of U, and let L OQ (x) = lim^oo Lt(x) for all x G R . By 
Lemma [3.121 and (a scaled version of) Proposition 12. 12l 

is not compactly supported ==>- U survives. 

It therefore suffices to show that is not compactly supported with positive probability. 
To do so, we will specify an algorithm that produces a (random) set £1 consisting of integer 
sites such that 

(i) L OQ (Q(x)) > for all ifSl; and 

(ii) f2 is infinite with positive probability. 

The set f2 will be the connected cluster containing the origin in a 3— dependent site perco- 
lation process with density > 1 — 6eo- 

Let us first give an overview of the algorithm. Recall that the grid T is defined to be Z+ 
when d = 2 and Z+ x {0} when d = 3. Initially all sites x G T are designated vacant (that 
is, Q = 0). Our algorithm relies on the comparison in Proposition 12.231 Starting from the 
origin, following the total order = x(l) -< x(2) -< ■ ■ ■ on T introduced in (j5.3D . we shall 
define stopping times n , random measures [n , i/j and suppression rates K* . Proposition 12.231 
allows us to couple U with another process U* , which, on any time interval between two 
successive stopping times, is a usual spatial epidemic process. The set Q will be determined 
by U*. Proposition 12.231 ensures that > , which will be used to ensure property (i). 
Depending on how U* behaves for t G [7"i_i,Tj], we may change the status of site x = x{i) 
from vacant to occupied, and add x to the set O. Roughly speaking, this will be done if 
and only if the spatial epidemic U% for t G [Tj_i,Tj] succeeds in (1) putting enough mass in 
adjacent cubes at time Tf, and (2) accumulating only a small amount of local time. On the 
event that the status of site x is changed to occupied, for each successor y G A{x) we will 
be able to extract a "nice" mass distribution fi y in such a way that if a spatial epidemic is 
initiated by [i y then it will have high probability of making the two events (1) and (2) occur; 
in other words, so that site y will also be added to O with high probability. By keeping this 
probability above the percolation threshold we will ensure that the random set consisting 
of all the occupied sites will be infinite with positive probability. 

We now introduce some notation. In addition to Oo, assume M, A" are as in Proposition 
15.41 so that (I5.1ip holds. In the algorithm, we will repeatedly use stopping rules r = 
t(Y;£;R) defined as follows: for a measure-valued process Yt G C([0, oo); M. c (M. d )) with 
local time LY , a threshold I > 0, and a region R C M. d , 

t(Y;£;R) :=mf{t : maxLf(x) > £ or Supp(Lf ) % R} A T. (5.12) 
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We will also repeatedly use the notation F 1 , . . . , F A as introduced in (|5.4j) to define the 
so-called "good events". For notational ease, for each % = 1, 2, . . ., associated with site x(i) 
in the above overview, we write XJ\ = U^ T . for t > (i.e., the process U* shifted and 
restricted to t > Tj_i), and 



G i = F\M; U\x{i)) n F 2 ( X ; U r ) n F 3 (U\ x{i)) n F 4 (e ; U 



). (5.13) 



The event G l will be what's call a "good" event. In plain language, ignoring the technical 
restriction F 4 , on such a good event, before time T, the spatial epidemic U l has not ac- 
cumulated local time density more than \ ' l<3~(a;(i))> an( i m the meanwhile, at time T, it 
spreads at least \Uq\ amount of mass in all the cubes Q(y) for y £ A(x(i)). 

Now we describe our algorithm in detail. In order to apply Proposition 12.231 we need to 
define four sequences: random measures suppression rate functions K* and stopping 

times Tj. The random measures /ij and will be defined through an auxiliary random 
measure sequence W{. The suppression rate functions K* will be deterministic functions as 
follows: Kq = 0, and for i > 1, 

i 

3=1 

Observe that for each i, K* is a summation of moving windows, and is bounded by (3n 
everywhere (recall that k is defined in (15. 9p ). 

We start with site x(l) = 0. The To, fio and vo are all deterministic: To = 0, = M 
and uq = 0. Let t\ = t(U*;XiQm( x (1)))- By Proposition 15.41 the good event G 1 occurs 
with probability > 1 — 3eo • Observe also that T\ > almost surely and t\ = T on G 1 . If the 
good event G 1 occurs, then we change the status of site to be occupied. Further define 



Wi 



E ff* 77T( X\ •^-(•HQW), if G 1 occurs; 
zeA{ x (i)) ^n-WWJ (5.15) 

k 0, otherwise. 



We now work with site y = x(i) for i > 2. We proceed according to whether the site y 
is an immediate offspring of some occupied site or not. 

Case I. Site y is an immediate offspring of some occupied site. Define 

(Hi-!, Vi-x) = n Q(y)), Wi -i(- n Q{yf)) . (5.16) 

Then is a measure supported by Q(y), of total mass \fi\ = fd(0), and (A(eq/A), fd(Q), fo( £ o 
admissible. Let Tj = r%-\ + t(U 1 ~ 1 ;x^Qj^(u))- By Proposition 15.41 (with an apparent spa- 
tial translation), the good event G % occurs with probability > 1 — 3eo- Observe also that 
Ti — Tj_i = T on G % . If the good event G l occurs, then we change the status of site y to 
occupied. Moreover, according to whether G l occurs or not, we define W{ as follows: 



Wi = < 



E n* iklTYs -C(' n Q( z ))> if occurs; 

otherwise. 
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where 

A(y) ={z£ A(y) : z ^ .A(u) for u which is occupied and -< y}. 

Case II. Site y is not an immediate offspring of any occupied site. Then we set (/ij-i, Vi—%) = 
(0,W{-i), Ti = n-i, and Wi = Wi-\. 

In either case at time Tj we proceed to site x(i + 1). 

It is easy to see that such defined [ii,Vi,K* and t% satisfy the conditions of Proposi- 
tion I2.23( and therefore the processes U and U* can be coupled such that 

Lf > I%* , for all t > 0. 

Now if we let f2 be the set of all occupied sites, then by the algorithm above, for any 
x = x(i) G f2, 

L£(Q(s)) > L^(Q(x)) - I%"_M X )) > °> 

and hence f2 satisfies the condition (i). 

We now show that fi is infinite with positive probability. Define a site percolation on T 
as follows: for each x G T, if x is occupied, then we let £(x) = 1 if both y G »4(:r) are 
occupied, and = otherwise; if x is vacant, then we let be a Bernoulli(l — 6eo) random 
variable that is independent of everything else. 

We know that the origin is occupied with positive probability. We claim that on the 
event that the origin is occupied, f2 contains the collection of sites reachable from the origin. 
We may assume that £(0) = 1 since otherwise we are done. But when is occupied, £(0) = 1 
implies that both y G .4.(0) are occupied. By induction the conclusion follows. 

It remains to show that the above defined site percolation is a 3-dependent site perco- 
lation with density at least (1 — Geo), i- e - 5 we need to show that for any n > 1 and any 
1 < i\ < . . . < i n such that — > 3, 

P(t(x(ij)) = for all j = 1, . . . ,n) < (6e ) n . 

Since when a site x is vacant, is a Bernoulli (1 — 6eo) random variable independent of 
everything else, we need only to show 

P (£(x(ij)) = for all j = 1, . . . , n \ all x(ij)'s are occupied) < (6eo) n . (5.18) 

Let us first consider the n = 1 case. When x := x{i\) is occupied, by construction, 
each y G A{x) is occupied with probability at least 1 — 3eo 5 hence the probability that both 
y G A(x) are occupied is at least 1 — 6eo- (|5.18p follows. 

In general, for each m > 0, we define Q m to be the cr-algebra generated by {C/ t * : 
< t < T m }. Then for each i > 1, the good event G l is measurable with respect to Qi, 
and hence the Bernoulli random variable £(x(i)) is measurable with respect to Qi where £ 
is the index of the second y G A(x(i)). Now since x(ij)'s are at least distance 3 from each 
other, if we let £j be the index of the second y G A(x(ij)), then 

£j < - 2, for all j < n. 

Hence by further conditioning on Gi„-2, (j5.18[> reduces to the n = 1 case and hence holds. □ 
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6 Proof of extinction when d = 2 or 3 



As the title suggests we shall assume d = 2 or 3 throughout this Section. 



6.1 Scaled process 

Proposition 6.1. Suppose Ut is such that for each ip £ C 2 

U t (i/>) = tfoW0 + ? f U s (A^j) ds + e I U.ty) ds-p fu s (L? ■ da+^yMtty), (6.1) 
1 Jo Jo Jo 

where Mt(ip) is a martingale with quadratic variation [M(ip)]t = Jq U s (^ 2 ) ds. There exist 
positive constants Eq and C such that if the initial condition Uq belongs to the class 

C := {/i satisfying Assumption \l.l[ Supp(^x) C Q(0), and = 2}, (6-2) 

and the positive parameters a, £,(3 and 7 satisfy As sumption 1 6. 2\ below, then 

P{U dies out) = 1. 

Assumption 6.2. 

P fa Js~\ , . fp /3 2 1 p\ „ 

e<——ji max £,-, — < e , and min [-^, -5-, -, - > C. 
2 • 3 a V 7 7 / \ £ £ 1 7 7/ 

We denote by p^ e ^ n the law of Ut satisfying (|6.ip with Uq = fx E C. Then we can 
rephrase the conclusion of Proposition 16.11 as 

where 



p a,e,P,<y ._ gup pa&PflQy survives ). ( 6 . 3 ) 



n a,e,/3,7 — gU p pa,e,/3,y 

When there is no confusion about the initial configuration /i, we omit /x and write p a > £ >P^ 
and sometimes just write P. Note that P a ' e >°'T denotes the law of a Dawson- Watanabe 
process without any local time killing, and _p Q > >°>T the law of driftless Dawson- Watanabe 
process. By (a scaled version of) Proposition 12.121 we see that when ft > 0, 

jja,e,P,l < {ja,e,0,7 ; (6.4) 

where C/ a,e '^' 7 has law P^' £, ^ n , [/~ a > £ >°>7 h a s J aw p^> £ >°>7 a nd the above notation means we 
can define versions of these processes on the same space with U^ Al < t/f* 6 ' ' 7 for all t > 
almost surely. Furthermore, by Lemma [2. 11 the laws p™' e ' ' 7 and P^' ' ' 7 are related to each 
other via the likelihood ratio 



, p a,e,0,7 

(U) 



dP^ 



a,0,0,7 



2 ft 



exp (-^ M t(!) " ^ / 1^1 ds ) . (6.5) 
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We introduce the following notation: 

V t = \U t \; r 3 = inf{t : Supp(Lf ) % Q 3 (0)}, 

and for any continuous real valued process Xt and any c 6 1, we let T C {X) be the hitting 
time: 

T C {X) = inf{i : X t = c}. 
Finally, define r to be the first time that Vt hits or 4 or that Ut exits Q3(0), i.e., 

r = T (y)AT 4 (y)Ar 3 . (6.6) 

Lemma 6.3. r < oo almost surely. 

Proposition 6.4. There exist constants Eq and C such that if the parameters a, e, (3 and 7 
satisfy Assumption \6.2l then 

sup P^(V T > 0) < := p c . (6.7) 

M 6C 2 • 3 a 

We will prov e these results in t he ne xt sub-Section. Proposition 16.41 is analogous to 
Lemma 2.3.1 in iMueller and Tribe! (|1994l ). Once we have the Proposi tion, we can prove 
Propo sition 16.11 by constructing a sub-critical branching process as in IMueller and Tribe 
()1994h . or more directly as follows. 

Proof of Proposition \6.1[ Suppose that the positive parameters a, e, (3 and 7 satisfy the as- 
sumption of Proposition 16.41 Let 



S up^ c P^(V T >0) ^ 



n 

Pc 



By the definition (I6.3P of p a ^ £ ^^ 7 we can find a fi G C such that 

P«> £ >^(U survives) > l±I!p^A7. 



Let Ut satisfy (16. lh with Uq = /x. For this U, at time r, on the event that V T > 0, U T is 
contained in Q3(0) with total mass no greater than 4. We can then decompose it into no 
more than 2 x 3 d parts: 

t 

U T = J2 U T> 1 < 2 x 
i=i 

each of whom has support contained in a unit cube, total mass at most 2, and satisfies 
Assumption 11.11 To see this last property, the domination in (|6.4p and absolute continuity 
in (16. 5p show that it suffices to prove that if U is the super-Brownian motion with law P^ ,0 ' ] 



it- > 

then for every fixed 5 > 0, U TS satis fies Assumptionll.ll a.s.. where 75 = (5 Vr) A5 _1 . Let us 
take d = 3. By Theorem III. 3. 4. in iPerkinsI (120021 ) we may fix oj outside a null set so that 
for some finite C = C(oj,8), 

sup U t (B(x, r)) < Cr 2 (log + - + 1) for all r > 0. (6.9) 

x,t>8 r 
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For any (t n ,x n ) -> (t,x), we want to show that J y qt n (y-x n )dU TS (y) -)• J y q t (y-x)dU TS (y). 
Fix an 7/ 6 (0, 1). The elementary bound ^(x) < c|ic| , (|6.9p and an integration by parts 
(left for the reader) lead to 

sup / q tn {v ~ x n ) 1+71 dU TS (y) < c sup / \y - Xnl" 1 ' 71 dU s (y) < 00 a.s. 

n J s>5,nJ 

Therefore {qt„{y — x n )} is uniformly integrable with respect to U Tg . We also have 

Qt»(y - x n ) -> qt(y - x) for U TS -&.a. y 

because U Tg ({x}) = by (16. 9p . The latter also shows that the limiting qt(y — x) is well- 
defined for U TS -&.&. y. Therefore we may take limits under the integral J qt n (y — x n )dU TS (y), 
and conclude that U Tg satisfies Assumption II .1| as required. A similar argument applies for 
d = 2. 

By the Markov property of the joint process (U t ,LY) (see Theorem 122(d)), Lemma^H 
and (a scaled version of) Lemma [2.171 

pofJnqj survives) < E (l(y T>0 ) ■ P^ Al {U survives)) . 

Here we are "throwing away" the killing due to . By Lemma [2.181 and translation invari- 
ance, the right-hand side is bounded above by 



a,e,/3,7 



E (l { v T >o) • Yl p ut A ^ U survives )) < p (Vr > 0) • E ^5> a ' s '^J < r lP 
Combining this with the previous inequality and (|6.8p we get 

hence p Q > £ '/ 3 >'>' = 0. □ 
6.2 Proof of Lemma 16.31 and Proposition 16.41 

In the arguments below, U% is a process satisfying (|6.ip with a fixed initial condition /x 6 C. 
The bounds in Lemmas 16. 6| - [rT9l below hold for all [J, EC, and hence will lead to the uniform 
bound in Proposition 16.41 

First we note that Vt = \UA satisfies the following SDE for some Brownian motion W: 



dV t = eV t dt - pU t (L? ) dt + V7 y/V t dW t . (6.10) 
By an integration by parts, 

P f U S {L U S ) ds = ^- I (L?(x)) 2 dx. 
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When t < T3, by Cauchy- Schwartz we get that 

^j\ L ^x)) 2 dx>^.^ d ^f^(x)dx^ 2 =p c ^j\ s ds^ . (6.11) 



We now prove Lemma [67 

Proof of Lemma \6.3l Suppose otherwise P(r = oo) > 0, in particular, P(t^ = oo) > 0. By 
(|6.10p and (|6.1ip . on the event {t% = oo}, 

V t < 2 + £ J V s ds-p c l3 (J V s ds^j + ypy J y/V s dW s , for all t > 0. (6.12) 
Define a sequence of stopping times {r{\ by rg = and for i > 1, 



n + i if V u +i < 2, 

inf{i > T{ + 1 : Vt = 2} otherwise. 



Claim 6.5. For all i, ri < oo almost surely 

Suppose for some i, T{ < rj+i = oo. Then Vt > 2 for all t > r, + 1. Therefore (16. 12j) shows 
that on {T3 = oo} the continuous martingale y^y J * x/Vg dW s approaches +oo as t — > oo, an 
event of probability zero. This proves the claim. 

For each i Proposition 12 . 1 2l allows us to bound Vt above on [fj, r, + 1] by a Feller diffusion 
with drift e and initial value 2 which does hit in the next one unit of time with probability 
q > 0. This shows P(V hits on [rj,rj + 1]|J>.) > q > and we therefore conclude that V 
will hit almost surely, again a contradiction to our supposition. □ 



a : 



Next we prove Proposition 16.41 Define a continuous random time change 
~0j^ V) V s ds] ^[0,T (y)]by 



a t = inf jr : ^ K da = t j , (6.13) 

and let Vt = V at . Then Vt satisfies 

V t = 2 + et-p U S (L U S ) ds + VrBt, for t < / F s ds, 
Jo Jo 

where 5 t = J Q at y/V s dW t for f < J Q To(V) V s ds and may be extended, if necessary, to a 
standard Brownian motion. If 

Y t = 2 + et - Pc (3t 2 + yfiBt. (6.14) 

then by ([BTTTl) 

V* < Y t for t < / V s ds, (6.15) 
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since the upper bound on t implies at < T3. 

We want to bound P(V^ > 0) where r is defined in (|6,6p . Using the comparison above, 
noting that by Lemma [6.31 t < 00 almost surely, we get that 

P(V T > 0) 

<P(r 3 < T (V),t 3 < T A {V)) + P(r = T 4 (V)) 

<P(r 3 < l/(4e)) + P(l/(4e) < r 3 < T (V)) + P(T 4 (U(- A 73)) < T (F)) 
<P(r 3 < l/(4e)) + P(Ti(V) < l/(8e) and T (V) > l/(4e)) (6.16) 
+ P(T (F) > l/(8e)) + P(T 4 (y) < T (Y)), 

where in the last line we used that 

PiT^V) > l/(8e),T (V) > l/(4e) and r 3 > l/(4e)) < P(T (Y) > l/(8e)). 

This holds because Vq = 2, and hence on the event on the left-hand side 



r 3 AT (V) 



V, ds > 



l/(8e) 



V s ds > 1/(8e), 



which implies av/sg) < l/(8e), and by (|6.15|) . for all t < l/(8e), Yt > Vt = V at > (since 
Po(X) > Proposition 16. 41 will be proved if we can show that all the four probabilities 

in (|6.16p are small. 

Lemma 6.6. There exists a constant C > such that 

P(r 3 < 1/(4.)) < C^exp (^) + 2^exp (^) . 

Proof. By the domination (16.41). it sufhces to s h ow th e lemma for which is then 

analogous to Lemma 2.1.9 in Mueller and Tribe (Il994h where the conclusion for the d = 1 
case is proved. We give here a slightly simpler proof for all d < 3. 



Following [Mueller and Tribe! (|1994l ) and using (|6.5p , we get that 

p^,0,7( T3 < 1/(4e)) 
<^ A0 ' 7 (W/WAT (V)} • ^ (^ M -(l) " ^ jH Vs ^ 

+ P Q ' £ '°> 7 (T E _ 1/2 (T0 < l/(4e)) 
<VP q '°'°' 7 (t 3 < l/(4e)) 



(6.17) 



/£ Q ,o,o, 7 ( exp ^-^M T3 (1 



V, ds 



7 Jo 

+ p^^(r £ - 1/2 (y) < 00). 



L{73<l/(4e)AT ' i /a (V)} 



(6.18) 
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A scale function (see, e.g., Proposition VII. 3. 2 and Exercise VII. 3. 20 in iRevuz and Yor 

:or V when (3 = is given by s(x) = 7(1 — exp(— 2ex / j)) / (2e) and so 



P a ' £ '°^(T e . 1/2 (V) < 00) 



s(2) - s(0) 1 - exp(-4e/ 7 ) 



< 



s ( £ -i/2) _ s (o) 1 - exp(-2^/ 7 ) 
4e/ 7 



(2v^/ 7 )exp(-2Vi/7) 



2 v / eexp(2 v / e/ 7 ). 



(6.19) 



We will use Theorem 1 of llscod ()l988l ) to bound P a A0,7( T3 < l/(4e)) < P a A0,7( T3 < qq). 
To do so, we make another scaling: let 



, for all tj> E C*{R d ), 



then by Lemma 12.261 Uf satisfies the assumptions of Theorem 1 in llscod (j 19881 ) , and 



U t (Ql(0))>0^U t lQl rj(0)\ >0. 



Hence by Theorem 1 in llscod (j 19881 ) and the fact that £7q E C 



P Q '°^(r 3 < 00) < ^ — 2 LL <Cn((2/3) ei )--, 



(6.20) 



3 2 

2 V a 



where e\ is a unit vector and u(x) is the unique positive (radial) solution of the singular 
elliptic boundary value problem 

Au(x) = u 2 (x), x E -8(0, 1); and u(x) — > 00 as |x| — > 1. 

To bound E a ^ (exp (^M T3 (1) - £ / Q r3 V s ds) ■ l{r 3 <i/(4 e )AT _ 1/2 (v)}) , we note that 
if Xt = Mj(l), then B s := X a t s \ may be extended to a standard Brownian motion, where 
a is the random time change defined in (|6.13p . Note that when T3 < l/(4e) A T £ -i/2(V), 
[X] T3 = £ 3 V s ds < l/(4e 3 / 2 ). Therefore 



E 



a, 0,0, 7 



exp 



X. 



2e 
2e 

7f~ 



<E [ exp ( _ j B 1/(4e 3/2 ) 



r 2 /-T 3 



7 Jo 



V.ds\-1 



L H<V(4e)AT -1/2 (V)} 



<^a,0,0, 7 ( exp ( 



2e 



7 

£ 



[ X ]^3 ) • 1 {[X] T3 <l/(4 E 3/2)} 



2 1 



7 4e 3 / 2 



: exp 



47 J ' 



(6.21) 
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where in the second to the last line we used optional sampling theore m is a stopping 

time f or the filtration {J 7 Qs }, see, e.g., Theorem 3.4.6 and its proof in lKaratzas and Shreve 
(|l99lh ). 

Now insert ([6TT9]) . (frT20|) . and ([OT]) into ([63?]) to complete the proof. □ 
Lemma 6.7. There exists a constant C > such that 

P(Ti(V) < l/(8e) andT (V) > l/(4e)) < C-. 

7 

Proof. Recall that satisfies (|6.10p . Applying Proposition 12. 12l again, on {Ti(V) < oo} we 
may define an Fx 1 (v)+t- adapted solution V of 



V t = l + e f Vsds + ^jl JV s dW' s , 
Jo Jo 



where W 1 is an J 7 ^ 1 (v r )+rBrownian motion and Vt > V^iy^t for all t > 0, almost surely on 
{Ti(F) < oo}. Therefore 

P(Ti(V) < 1/(8e) and T (F) > l/(4e)) < P(T (F) > 1/(8e)). 

By Exercise II. 5. 3. in IPerkind (|2002h the last term equals 

/ -2e \ 2e 
1 " ex P -7; TTT^T < 



7 (l_ e -i/8)y - 7 (i_ e -l/8)" 

□ 

Lemma 6.8. There exist constants C\ y Ci > such that for all (3e~ 2 > 20000 and < 
e< 1/4, 

/3 2 



P(T (y) > l/(8e)) < Ci exp -C7 2 ^- 

V e 7 

Proof. Assume /3,e are as above. Recall that Y is defined in f|6. 14|i . 

P(T (Y) > l/(fe)) < P(yi /(8e) > 0) 

1 /Pc „-2 - 1 



<Fk>^(r 2 - ioooo) 



- ^ X " 128 ^7 1 V 
The result follows. □ 
Lemma 6.9. There exists C > suc/i that if e < min (1/2, p c /3), then 

P(T 4 (Y) < T (Y)) < Cexp + exp ("2p/ 
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Proof. Recall that Y t satisfies 

Y t = 2 + et-p c f3t 2 + ^B t . 

Hence if we define Y t by 

Y t = 2 + et + ^B t , 

then Yt < Y%. Note also that 

P(T 4 (Y) < T (Y)) < P(T 3 {Y) < 1) + P(T 3 (Y) > 1,T 4 (Y) < T (Y)). 
We first estimate P(T 3 (Y) < 1): 

P(T 3 (Y) < 1) = P(max(et + ^B*) > 1) 

1-e 



< P max B t > — — (a 99^ 

- \t<i l - v^J (6 } 



< Cexp (- — 

V 87, 

provided that e < 1/2, where C > is some constant independent of e and 7. 
We now work with P(T 3 (Y) > 1,T 4 (Y) < T (Y)). Define 

Y t = Y 1 -p c pt + y/7B t , 

where B t = B t +i — B\. If e < p c /3 and t > 1, then 

Y t = Y 1 + e(t - 1) - p c /3(i 2 - 1) + ^B t -i 
= Y 1 + (t-l)[e- p c (3{t + 1)] + yfiBt-! 
< Y 1 + (t - l)[p c /3 - 2p c /3] + VTB t -i 
= F t _i. 

Furthermore, since F t < F t , on the event {T 3 {Y) > 1}, Yi < 3 and T 4 (Y) > 1. Therefore 

P{T 3 {Y) > 1, T 4 (F) < T (Y)) < P (3 - p c /3t + ^7^ hits 4 before 0) . 
The latter probability can be explicitly calculated using scale functions: if we let 

> c /3 



-x 



1 



then 



Vcf> ^bu-. „\ »(3)-»(0) «p(2^-S 



F|3-^. t+ «M ts 4W™ 0/ s(4) _ s(Q) ^ 

<exp \-2p c - 
V 7 



□ 
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Proof of Proposition \6.4\ The hypotheses of the above 4 Lemmas are satisfied under As- 
sumption 16.21 for small enough £o and large enough C. The bounds obtained in all four 
Lemmas can also be made as small as we like, again by taking Eq small enough and C large 
enough. By inserting these bounds into (|6.16p . we obtain Proposition 16.41 □ 

6.3 Proof of extinction for the original equation 

By Proposition 14.31 and Proposition 16. 1\ in order to show extinction for Xt defined by the 
original equation (jl.ip . it suffices to show that when 6 > is sufficiently small, there 
exists a scaling as in Lemma 12.261 such that the parameters in the scaled equation satisfy 
Assumption 16.21 This is the content of the next lemma. 

Lemma 6.10. For any fixed constants < Eq < C , for all 6 > sufficiently small, there 
exist a scaling of X, as in Lemma \2.2b\ with K = 0, such that the parameters in the scaled 
equation satisfy A s sumption 1 6. 21 

Proof. By Lemma [2.261 we want to find positive constants a, b and c such that 

a 2 b d 

a = ab 2 , e = a6, (3 = , and 7 = ac (6.23) 

c 

satisfy Assumption 16.21 We will only look at power functions, i.e., 

a = e x ,b = 6 y , and c = 6 z , 

and show that for appropriate (real) choices of x, y and z, Assumption ^. 2l is satisfied provided 
that 6 is sufficiently small. We have that 

a = 6 x+2 v, e = 9 1+x , p = e 2x+d y~ z , and 7 = 9 X+Z . 

Looking back at the conditions in Assumption 16.21 we see that it is sufficient that 

' 1 + x > 2x + dy — z, 
1 + x > 0, x + 2y - (x + z) > 0, (1 + x)/2 - {x + z) > 
(2x + dy- z)- 2(1 + x) < 0, 2(2x + dy - z) - 3(1 + x) - (x + z) < 0, 

, x + z < 0, 2x + dy — z — (x + z) < 0, 

i.e., 

' 1 + z > x + dy, 
x > -1, 2y > z, 1 - x > 2z 

< 

dy - z < 2, 2dy < 3z + 3 
, x + z < 0, x + dy < 2z. 

There is an abundance of choices, for example, x = —3/4 and y = z = 1/2 will do. □ 
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7 A strong form of local extinction 



Theorem 11.31 is immediate from Proposition 14.41 and the following result. 

Theorem 7.1. Assume that d = 2 or d = 3. If the initial mass distribution /i satisfies 
Assumption li.il then for any value of 6 the epidemic Xt (the solution to the martingale 
problem (1.1)) dies out locally, that is, with probability one, for every compact subset K C R rf , 

Xt(K) = for large enough t. (7.1) 



The remainder of this section will be devoted to the proof of this theorem. Observe at 
the outset that it suffices to show that the property (|7.ip holds when K is a ball of radius 
g = g{9) > centered at a point with rational coordinates, because any compact K is 
covered by finitely many such balls. Moreover, it suffices to consider only balls centered 
at the origin, because the initial mass distribution ^ can always be re-centered. Thus, our 
objective is to prove that the epidemic dies out in K = B g (0). 



7.1 Re-infection at large times 



The proof of Theorem 17.11 will have three parts: First, we will show that (|7,ip could fail 
only if the ball B e (0) were re-infected from outside the ball i?3 e (0) at indefinitely large 
times. Second, we will show (in Section 17.21 below) that boundedness of of EL oo (B3 S (0)) : by 
Proposition I4.1[ implies that the mean mass flux through the sphere of radius 2g is finite. 
Finally, we will show (in Section 17. 3ft the finite total mean mass flux through the sphere of 
radius 2g will imply that reinfection of B g (0) from outside ^3^(0) at arbitrarily large times 
cannot occur. 

To give precise meaning to the notions of "re-infection from outside" and "mass flux 
through a boundary" we must bring in the historical process Ht associated with the spa- 
tial epidemic Xt. (For a rigorous development of the b asic t heory, for Dawson- Watanabe 
processes without inter action, see iDawson and PerkinsI (|199ll ). for interac tive proc e sses i n- 
cluding our setting see Perkins (j 19951 ). and for an overview of both see PerkinsI ( 2002 )). 
Recall that for each time t the state Ht is a random measure on the space of continuous 
paths C([0,t],R d ) that projects to X t via the time— t evaluation mapping. As in the above 
references, for w G C : = C([0, oo),M d ) we set w t (-) = w(- A t) and identify C([0,t],R d ) with 
{w E C([0,oo),M d ) : w = w 1 }. 

Theorem 5.11(a) of PerkinsI ( 1995 ) gives a version of Dawson's Girsanov Theorem for 
historical proces ses. It i s then easy to adapt the proof of Theorem 12.21 to see that The- 
orem 5.11(a) of PerkinsI ( 1995 ) will apply with the drift function g there equal to 6 — 
Lf(w s ). This gives a solution Ht to a well-posed historical martingale problem so that 
Xt((p) = f ip(wt)Ht(dw) is the unique solution to (jl.ip . It also shows that the law of H is 
absolutely continuous to the law of the historical process associated with super-Brownian 
motion on the filtration up to time t, for each t > 0. 

For a fixed g > let 



rj{w) = rj e {w) = inf{t > : \w%\ > 3^} 
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be the exit time of the path w from the interior of B3 e (0). At time t color the path (w s ) s <t 
red if r] < t and otherwise color it blue. This gives a decomposition 

H t (.) = H t R (.) + Hf (•) := H t {- n {q < t}) + H t (- n {r? > t}). (7.2) 

Projecting via the time— t evaluation we obtain the decomposition 

X t (.) = Af (•) + X t B (.) := H«(wt € •) + H t B (wt € •)• 

Proposition 7.2. For eac/i value 9 € R i/iere exists g = g(8) > suc/i i/iai /or any initial 
mass distribution \i satisfying Assumption \l.ll the process H B in the red/blue decomposition 
(|7.2p will die out with probability one. 



Proof. Arguing as in Proposition IV. 1.4 of Perkins! (j2002l ). but using historical processes, 



one can construct our historical epidemic process H and the historical process H for a drift-# 
super-Brownian motion, X, on a common probability space so that Hq = Hq and Ht < Ht 

Q 

for all t > 0. We decompose H = H + H as in (|7.2p . thus inducing a corresponding 

Q Q 

decomposition, X = X + X . Then X will be the drift-0 superprocess associated with 
Brownian motion killed when it exits the interior of B^ g (0). Therefore if Q x is Wiener 
measure starting at x and rj e is also the corresponding exit time for the Brownian path, then 
for t > 0, 

E(\xf\) = e et J Q x (t < ri e )dn{x) < e e %\Q (t < r, Q ). (7.3) 

(A careful proof of this could use the appropriate version of Proposition 17.4( c) below with 
r\ e in place of and H in place of H.) For g > sufficiently small (how small will depend 
on 9) this expectation decays exponentially with t, by elementary estimates on the transition 
kernel for killed Brownian motion. (In particular, 3^ > must be small enough that the 
first eigenvalue of —A/2 with Dirichlet boundary conditions on dB^ (O) is strictly greater 
than 9.) 

Q Q 

It remains to show that the exponential decay of E\X t \ implies that X dies out almost 
surely. Let Z± denote a Feller branching process with drift 9. For n G N, the fact that the 
total mass process \X | is dominated by the total mass process without killing on dB^ Q {{)) 
implies 



P{\X B n+1 \ > 0) < E{P-b (|xf | > 0)) < E(P { -b >0))=E 



Q 

, ) -^\x n \ 

1 — exp 



< C{0)E{\X*\) t (7.4) 



where Exercise II. 5. 3 of PerkinsI ( 2002 ) is used in the next to last line. The exponential 



decay in the mean on the right-hand side now shows that H B , and hence the smaller H B , 
dies out a.s. by a Borel-Cantelli argument. □ 

For future reference we state a time shifted version of the above. Let T > 0, define 

a T = inf{t > T : \w t \ > 3g}, (7.5) 

and for t >T set 

H t B ' T (.) = H t (.n{a T >t}). 
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Proposition 7.3. For fj,, 9 and g(9) as in Proposition 1 7. 2[ the process H t ' will die out 
with probability one. 

Proof. One proceeds just as above but conditional on the past up to T, H t , t > T will be the 
historical process associated with a drift-# super-Brownian motion starting at Ht = Ht- □ 

Assume for the remainder of the proof that g = g(9) > is small enough that the 
conclusions of Propositions 17.21 and 17.31 hold. Then for any fixed T > 0, all mass in the 
spatial epidemic Xf will eventually be descended from the mass in Xt outside of £?3 e (0). 
This obviously implies that if local extinction 07.1 P fails for K = B g (0) then the ball B g (0) 
must be re-infected by mass from outside B^ g (0) at arbitrarily large times. 



7.2 Finite mass flux 



We will control the re-infections of B g (0) from outside i^^O) by bounding the total "mass 
flux" (to be made precise below) through dB2 g (0). For any continuous path w in M. d define 
vq < T\ < v\ < ■ ■ ■ to be the successive times of passage between the spheres dB-^ e {Q) and 
dB2g(0) (that is, vq is the first hitting time of B^ g (0), T\ the first hitting time of dB2 S (0) after 
uq, and so on). Now for each k = 1,2,... define H^ to be an associated historical process in 
which historical mass frozen at time r/% is collected as occurs for Tk < t. For gen eral super- 
processes these are the historical random measures constructed by Dvnkinl ( 199ll ) (Theorem 



1.5) u sing log Laplace equations. We will follow Theorem 2.23 and Remark 2.25 of Perkins! 



(|1995l ) which gives a recipe for their construction and associated stochastic analysis, using 



historical stochastic calculus, and does so in a more general interactive framework which 
includes our spatial epidemic processes. 

(R d ) denotes the space of bounded continuous functions on IR^ with bounded contin- 
uous partials of order 2 or less. 

Proposition 7.4. For each k E N, there is a non- decreasing continuous Ai(C)-valued pro- 
cess, H k , and hence an associated random measure on [0, oo) x C (also denoted by H k ), 
satisfying Hq = and the following properties: 

(a) w = w Tk , Tk(w) = t, and so wt S dB2 g (0), for H k — a.a. (t,w) a.s. 

(b) If ip is a bounded measurable function on C, then with probability 1 for all t > 0, 

^)l {t>Tk) H t {dw) = jf j VK fc )l(s>. fc) dM H {s,w) 

+ f I ^ wTk ^>r k )\° ~ L*(w s )]H s (dw)ds + Hj°ty), 
where M H is the orthogonal martingale measure associated with H. 

(c) If Xt(-) = f J l(yj T e .^H k (ds,dw) and ip € C^(M. d ), then with probability 1 and for all 
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t > o, 



ip(w s )l( s>Tk ) dM H (s,w) 
Aip 



+ 



Hs>T k ) 



(w s ) + <p(ws)(6 - L? (w 8 )) H s (dw)ds 



L 2 



(d) For any fixed t > and bounded Borel ip : C — > M ; if 

oo „ 

A n (t,i/>) = J^l(a-»<i) y ^(^ Tfc )l((i-i 



)2- n <T k <i2~ n )Hi2 



i/jen A n (t,i/j) — > H^(i/j) in probability as n — > oo. // A n a/so denotes the measure on 
[0, oo) x C associated with A n (t+,t/j), there is a subsequence {rij} so that ^4 n |[0,r]xC ~~ ^ 
in .M([0,T] x C) for all T > a.s. 

Proof. The above result is implicit in Remark 2.25 in Perkins! ( 1995 ) and carried out for the 
total mass in Theorem 2.23 of the same reference. We will sketch how the latter construction 
is easily extended to the measure- valued process H k . 

Let i jj > be a bounded Borel function on C and in the setting of Theorem 2.23 in 



Perkins! (|1995h . set 



C(t,u,w) = ^(ty Tfc )l( t>rfc ( to )). 



(7.6) 



The above setting includes our historical epidemic process with the function g on page 9 of 
this reference equal to — Lf(uj,w s ) and the integrator Z° on page 12 given by 

dZ°(s,w) = dM H \s,w) + 6H s {dw)ds - L* (w s )H s (dw)ds. 

Therefore for ip fixed, Theorem 2.23 in Perkins ( 19951 ) implies (b) and the first conclusion 
in (d) for some non-decreasing left-continuous process H^(ip) satisfying Hq(i/j) = 0. To 
derive (c) from (b) (with ip{w) = ip(w Tk )) we need to show that 



(<p(wt) - ip(w Tk ))l (t>Tk) H t (dw) 



+ 



{<p{w s ) - <p(w Tk ))l {s>Tk) dZ°(s,w) 
^f(w s )l (s>Tk) H s (dw)ds, 



and this follows easily from the historical stochastic calculus in Chapter 2 of lPerkins! (Il995h . 

Consider next the continu ity of H^{ijj) in t for ip > bounded and Borel. By (IV. 48) of 
Dellacherie and Meyer ( 19821 ). it suffices to show that if T n \, T are bounded (J^-stopping 



times, then 



lirn £«W-4W) = 0. 



Arguing as in (2.44) of iPerkinsI (j 19951 ) this reduces to showing 

lim E{H Tn (T< Tk <T n )) = 0, 



(7.7) 
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and 

lim E(H s (T k = T)) = for each s > 0. (7.8) 



n— >oo 



We consider only (|7,7p as the proof of (|7.8j) will then be clear. Using the weak continuity 
of H one easily sees that 

lim sup E(H Tn (T <r k < T n )) < E(H T (\w T \ = 2g,r k < T)) < E{H T (\w T \ = 2< ? )1 (0<T) ), 

n—toc 

(7.9) 



where we used Tk(w) > 0. Theorem III. 5.1 of iPerkinsI (|2002l ) and our absolute continuity 



of X with respect to super-Brownian motion show that 

P{H t {\w t \ = 2g)>0 for some t > 0) < P(X t (dB 2g (0)) > for some t > 0) = 0. (7.10) 



This implies the right-hand side of (17.9j) is zero and so (17. 7j) is proved, thus giving the 
continuity of H k (tp) for each tp as above. 

Next we construct H k as a measure- valued process. Choose a countable determining 
class T> of bounded continuous functions on C containing the constant 1. For each ip £ D 
there is a subsequence {n^} so that 

sup\A n -(t,^) - H]?(if>)\ = for all T > a.s. (7.11) 

This holds by the first part of (d), monotonicity in t and the a.s. continuity of the limit. 
By diagonalization we assume the same subsequence works for all if) S T>. It is then easy to 
check that A n . |[o.t]xC ~^ -^ fe |[o,r]xC as finite measures o n [0, T] x C for all T > 0. Formally 
we may use Jakubowski's theorem (Theorem II. 4.1 of Perkins! ( 2002 )) and the fact that 



the required compact containment condition follows easily from the modulus of cont inuity 
for the historical paths of super-Brownian motion (Theorem III. 1.3 of Perkins ( 20021 )) and 



the usual absolute continuity argument. Implicit in the above notation is the fact that the 
limiting random measure H k is related to the processes H k {ip) constructed earlier by 

J J 4)(w)H k (ds,dw) = H?(ip) for all t > a.s. 

This gives the existence of the required process H k satisfying properties (b)-(d). 

We have Tk(w) < t, w = w(- A r^), and so Wt € dB2 S (0) for A n (dt, dw)-&.s., and taking 
weak limits in n we obtain (a) except with T k < t H k -&.s. To see that T k = t H k -&.s., it 
suffices to fix t > e > and show H k ((t — £,t] x {t/u < t — e}) = a.s. This is easily derived 
from (b) with ip = l( Tfc<t _ e ) and a bit of historical stochastic calculus. □ 

We may repeat the above construction with minor changes for the stopping times v k in 
place of Tfc and so obtain continuous non-decreasing .A/f (C)-valued processes {H k : k E N} 
and their projections {X k : k 6 N} which are .M (Revalued processes supported on dB3 (0). 
We identify X k and X k with the corresponding random measure on [0, oo) x M. d . 

For future reference we state a truncated version of Proposition 17.4( c) . For t > T > 
define 



X k ' T (-) = j T J l ( ^ fc6 .) l {T < Tk) H k (ds,dw). 
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Proposition 7.5. IfT > and ip G Cf (M d ) ; f/iera urat/i probability 1 for all t > T: 



(p{w t )l {T < Tk<t) H t (dw) 
t 



l(T<T k <sMw s )dM H (s, w) 
Aip 



+ J I 1 (T<r k <s) 



+ <p(0 - If) (w s )H s (dw)ds + (X t fc - 



L 2 



(7.12) 



^{w t )\ T <Tk,v k <t)Ht{dw) 

1 (T<T k ,u k <sMw s )dM H (s,w) 



(7.13) 



+ 



L (T<r fc ,i. fe < S ) 



+ <p(6-L?) {w s )H s (dw)ds + X^Xv) 



k,T, 



Proof. For (|7.12p start with Proposition 17.4( b) with tf)(w Tk ) = (p(w Tk )l^x< Tk ) and then 
proceed as in the derivation of (c) above. The fact that Tfe = s for H k -&.& (s, w) is used to 
get the form of the final term. The proof of (|7,13p is similar. □ 

The total flux measure on [0, oo) x dB2 g (0) is X T = Y^kLiX k and similarly we define 
X v = X^fcLi X k on [0, oo) x £?3 e (0). At present these measures may be infinite. 

As was already noted, our plan is to control the re-infections of B s (0) from outside 
B3g(0) by bounding the total flux, \X T \, through dB2 g (0). We next bound this flux in L 1 
as a consequence of Proposition 14.11 and Proposition 17.4( c) above. 

Color a path w yellow at time t iff Tfc < t < vj., for some k > 1, that is, iff at time t w is 
engaged in an excursion from dB2 g (0) to dB^ g (0). Let Hj be the restriction of Ht to the 

yellow paths at time t, that is, Hj {A) = f 1a(w) Y1T=i ^-(r k <t<u k ) H t (dw), and let Xj be 
the corresponding time-t projection. 

Proposition 7.6. E(\X T \) < oo and E(\X"\) < oo. 

Proof. We only prove the first conclusion as the proof of the second is similar. 

By differencing the decompositions in Proposition 17.4( c) for times and we see that 

for ip G C?(R d ), 



v{ w t)\ Tk <t<u k ) dH t {w) 



^{Ws)^{r k <s<u k ) dM H (s,w) 



(7.14) 



+ / n / l( Tk <s< Vk )(^(w s ) + <p(w s )(e-Lf(w s ))^H s (d W )ds + 
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Let < ifo < 1 be as above with support in the interior of B3 g (0) and so that cpo = 1 on 
£?2g(0). Then X^(ifo) = and Xf((fo) = \X^\ for all k, t. Take expectations in the above 
with (p = (fo and then sum over k to conclude that 

E(Xj{^)) = e(J J ^pK) + Mws)[0 - L?(w.)] Hj(dw)ds) + E(X T ([0, t]xC). 

Rearrange the above and use Xj < X t and then CQ) to see that 

E(X T ([0,t] x 0) < E(X t (<p )) + E(j\ s (^.+<p (LX + e-) > ) ds) 

< ufa) + E({L? ,\A<p \ + Vo0 + )). 

The right-hand side remains bounded as t — > oo by Proposition 14, II and so the result follows. 

□ 

7.3 Local extinction 

Recall that rj(w) = T] s (w) is the exit time of w from the interior of ^^(O). For any path 
w € C, if rj < t and \wt\ < 2g, then for some k > 1, < t < v^. That is, if you exit from 
the interior of B^ s (0) before time t and at time t are back in the interior of B2 g (0), then 
t must fall in one of the excursions from dB2 P (0) to dB^ p (0). Therefore if tp\ € C^°(M. d ), 
takes values in [0, 1], has support in Br 3 / 2 ) e (0), and ip\ = 1 on B g (0), and T > 0, then for 
all t>T, 

/OO p 
<Pi(w t )l(t <r ,)Ht{dw) +^ / Lpi{w t )l( Tk< t<v k )H t {dw) 

k=l 

oo „ 

= xf((p!) + ^2 / n(m)i(,T k <T<t<v k ) H t(dw) 

k=l J 

oo „ 

+ ¥'l( U; t) 1 (T<r fc <t<^)- H t(^) 
jfc=l ^ 

:=X t B (^) + lf' T (^) + xP(^). 

We have decomposed X y according to whether or not the fcth return to i?2 e (0) occurs before 
time T or after it. 

We have already shown (Proposition IT.2[) that X B dies out a.s. Recall the ot defined 
in (|7.5p . Clearly < T < t < implies ot > t for Hf&.a,. w for all t > T a.s. (recall 
(|7.1U|) ). and so by Proposition 17.31 

X t ' (tpi) < / ip\(wt)H t ' (diu) = for large t a.s. for each T > 0. 

Therefore to complete the proof of Theorem 17.11 it suffices to show the following: 
Proposition 7.7. lim^-^ P(X^ ,T (ip 1 ) > for some t > T) = 0. 
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To prove this result we first recall a standard method used to compute hitting probabil- 
ities for a super-Brownian motion Xt with drift 9. For A > let U(t,x) = U x (t,x) be the 
unique non-negative solution of 



dU _ A 
~dt ~ ~2~ 



Ut + eut-ufft + xtp!, u = o, 



(7.15) 



which is bounded on [0, T] x M. d for all T (e.g., see Theorem 1 1.5.11(b ) in Perkins! ( 20021 )). 
The duality for superprocesses (e.g., see Theorem 11.5.11(c) in iPerkinsI (|2002l )) implies that 
for all initial measures v, 

£„(exp(-A J X s (^)ds)) = eM-"(U t x )). 



It follows that U x (t, x) increases as A, t — > 00 to a Borel function U°°(x) > satisfying 



Pv(X s (<pi) > for some s > 0) = 1 - exp(-v(U°°)). 



(7.16) 



Next use the fact that X propagates locally at a finite rate (see Theorem III. 1.3 of Perkins! 
(Hoo2)) and dies out in small time with high probability if \Xq\ is small (recall (|7.4p ). to see 
that for a small, 

sup P a s x (X s (fi) > for some s > 0) < -. 

|x|>2£ 2 

It therefore follows from (|7.16p that 



sup U°°(x) = C p < 00. 

\x\>2q 



(7.17) 



Proof of Proposition 1 7. 71 Fix T > 0. By differencing the decompositions in Proposition 17.51 



we have for ip £ C^(R d ), with probability 1 for all t>T, 



ip(w t )l(T< Tk<t <v k )Ht{dw) 



(7.18) 



+ 



L(T<r fc <s 



l (T<T k <s<v k ) 



< Uk )(f(w s )dM H (s,w 



^f + <p{0- L?j\ (w s )H s (d W )ds + [X t k - X$\(<p) - X^ T {<p). 



Fix u > T. Arguing as in Proposition II. 5. 7 of IPerkinsI d2002h it is easy to extend (|7.18p 
to time-dependent test functions on [ 0, u] x M. d , including V(t, x) = U x (u — t,x) (see also 
Theorem 11.5.11(b) of IPerkinsI (2002J) for the regularity of the above V). One gets an 
additional term involving and so with the above choice of V, the equation (|7.15p shows 
that the function in the square brackets in the second integral in (|7.18p becomes 



dV AV y 



V 2 
* s 

2 



Xtpi 



L*V S . 



57 



Therefore for T < t < u, 

V t (w t ) 1 (T< Tfc <t<v k )H t (dw) 
J^V s {w s )l {T < Th<s < Uk) dM H (s,w) + J^ J : 

{w s )H s (dw)ds + J J V s (x)[X k - X k ' T }(ds,dx). 



1 (T<r k <s<v k ) 



V 2 



Rearrange the above and sum over k (using Proposition [7]6]) to see that if X U,T = X^fcLi X+' T , 
then for T < t < u, 

xF(Vt)+ f \X^ T { Vl )ds (7.19) 
Jt 

rt ,t/2 , r t r 



(ds,dx)-X u > 1 (ds,dx)}, 



where M t ' (V) is a continuous martingale, starting at at time T, and satisfying 
(M Y ' T (V))t = Jq xJ' T (Vg )ds. Using Proposition 17.61 we see the last term is continuous 
in t and it then follows easily that each of the terms in (I7.19P is continuous. Now apply 

—X t ' (U£-t) — A J T X s ' ((pi)ds I and take expectations at t = u, where 
V u = Uq = and note that X T ' = to deduce that 

E(l-exp(-X f X s y ' T ((^i)ds)) 

= E(J\xp(-X? T (Utt) -\f T X^ T (^)ds) 

x J U^_ t (x)[X T (dt, dx) - X u > T (dt, dx) - Lf \x)XJ' 7 \dx)dt] 

Let u, A f oo and drop the last two negative terms to show that 

X^ T ( Vl )ds > o) < E^J J U°°(x)X T (dt,dxf) < C g E(X T ([T,oo) x R d ). (7.20) 

The bound ([7.17P on U°° for |x| > 2p is used in the last inequality. If we sum (|7.18p over k 
we may argue as in the analysis of (17.19P to see that X^ ,T ((pi) is continuous in t. This and 
the fact that the upper bound in (I7.20P approaches zero as T — > oo by Proposition 17.61 imply 
the required result. □ 
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